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CHAPTER I 

VV " ^ 

INIRODUGTION-' 

This manual describes the necessary preliminary electric circuit imodei ing ' \ 
and methods of computer analysis using the NASAP digital computer program to 
aid in the design of control circuits used in aerospace systems. The scope of 
design treated is limited to single -input single-output linear time invariant 
control systems that can be described by rational transfer functions. 

Two different approaches can be followed in the design of such linear 
time invariant feedback control systems. In one case the designer seeks 
compensators for a given plant to satisfy the over-all system requirements. 

This approach generally involves cut and try and the well known conventional 
design techniques such as the Bode and Nyquist plots and use of the root 
locus. In the other case, the designer starts by obtaining the over-all trans- 
fer function from the given plant and the specifications. He is then in a 
position to determine the required compensators. The second of these two 
approaches seems more amenable to computer-aided design of control circuits 
using the NASAP program. 

HASAP is an acronym for Network Analysis for Systems Applications ProgramV' 
NASAP has been developed and is being maintained by the Automated Techniques 
Branch, NASA Electronics Research Center, Cambridge, Massachusetts, The NASAP 
'program is based on Mason r s signal-flow graph [MAd] as extended by Happ for 
the closed signal-flow graph [HA 1]. The program is based upon symbol oriented 
techniques which with proper tagging and loop evaluation permit both the trans- 
fer function and the sensitivity to be made available. 

This chapter provides brief discussions on the basic- options presently 
available in the NASAP program and the basic procedures used in computer aided 
design of control circuits. 




r^ ne^MSAP Program 
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MSAP is a linear electrical circuit analysis program which computes a . 
specified transfer function in terms of the complex frequency variable, s. 

The basic coding rules and description of the program are contained in the book= ~ 
let "Coding Instructions for MSAP 69 / 1 " by Gaertner Research Incorporated [GA l] 
written under MSA contract MS 12-663. A number of special options are,- -intro- - 
duced and described -later in this manual. The input for MSAP -isf by • punched card 
and the output is\Syprinter tabulation of data and printer graphics. It is 
worth noting that the input format of MSAP is simple and that use of this 
program requires little knowledge of circuit theory or computer programming. 

The present version of MSAP can handle linear circuits which consist of 
constant- value passive elements, and independent or dependent current and 
voltage sources. The dependent sources must be linearly related to a voltage 
or current in another part of the circuit. Nonlinear functional relationships 
(dependenciesj^and time -varying parameters cannot be handled. 

Embodied in MSAP is the ability to give both a mathematical formulation 
and a numerical tabulation of the output results with some printer graphics. 

A brief summary of the options presently available in the program available on the 
:_ar<®;Sp at the Moore School, follows: 

1 . OUTPUT - The transfer function which is specified by the user is 
printed out as the ratio of two polynomials in the complex fre- 
quency variable, s. The poles and zeros of the transfer func- 
tion are evaluated. 

. 2r. EREQ - A' Bode plot of the transfer function is printed in 
tabular and graphical form. 

3. TIME - The ’impulse^: response of the network is printed in tabular 
and graphical form. To facilitate control circuit design we 
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have also made available the step and ramp responses. 

4. SENSITIVITY - The sensitivity of the transfer function is com- 
puted with respect to a ^designated element in the network. Fur- 
thermore the program can print out in numerical form: the sensi- 
tivity of the real part of the transfer function, of the ima- - 
ginary part of the transfer function, of the magnitude of the 
transfer function, and of the phase of the transfer function 
to changes in specified circuit parameters. 

* Details of these options are documented in Section ID. 

Pfc^edures Used in. Computer Aided Design 

The major steps in computer aided circuit design are: 

a. Choose the electrical circuit model for the control system. 

b. Calculate element values from design equations. 

c. Analyze circuit using computer. 

d. Change element values if computer analysis results in 
discrepancies between the electrical circuit model response 
and desired control circuit responses 

Before a control system can be analyzed using NASAP and before the 
necessary compensation can be determined, the plant* s dynamic characteris- 
tics must be first simulated by an electric network which has an equiva- 
lent dynamic characteristic . The choice of circuit models is based on 
the specified response for the control circuit. Details of circuit modeling 
- are given in chapter III . Computer analysis is amply demonstrated in 
chapter IV. 



'• Jk - - REVIEW OF CODING PROCEDURES FOR 3M.SAP 
Introduction 

NASAP is a computer-aided electrical circuit analysis program which can 
be used by engineers without knowledge of computer programming. Using 
MSAP to analyze a circuit it is possible to calculate a transfer function, 
the sensitivity of the transfer function, the frequency response, and the 
impulse response. The main steps for preparing NASAP computer instructions 
are: 

a) Obtain an electrical circuit model for the control circuit 

b) Preparation of the circuit diagram in a form from which all 

the information required by the computer can be readily extracted. 

c) The preparation of the computer instructions themselves. 

The first step will be discussed in chapter r fli. The second .step -has been de: 
cribed elsewhere [GA l] and mil be illustrated by many examples throughout 
this manual . The third step will be reviewed in this chapter for the conve- 
nience of the reader. For additonal information, see [GA 1]'. The coding rul< 
will be presented by an example taken virtually unchanged from [MO l]. 

Example • RLC 7 Filter' . • _ 

f 

We shall analyze the RLC filter shown in Fig. 1.1(a) to find its open 
circuit transfer voltage function. 

Numbering the Nodes and Elements 

The first step in preparing the circuit for MSAP analysis is to number 
the circuit nodes and elements. The nodes must be numbered sequentially, 
starting with 1, without skipping any numbers. 

The second step is numbering all circuit elements (resistors, capacitors, 
inductors, and current and voltage sources) consecutively within each category 
i.e., Rl, R2, . • ., II, L2, . . Cl, C2, etc. No two elements can have 
the same number designation regardless of their component values. One 

- 4 - 
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N AS AP PROBLEM (RLC FILTER) 


m j 

VI 1 2 loO 
LI 3 4 1 .OH 

Rl 2 3 0*66667 

R 2 4 5 0.55556 

Cl 5 l 0.36F 
R 3 4 1 0.8333 

^output 

VR3/VV1 • 

FREQ -1*2 1.0 0 
TIME 5.0 

^execution 


COLUMNS 


NETWORK 
Y ELEMENT 
TABLE 


05 


OUTPUT 

SPECIFICATIOWS 


CONTROL 

STATEMENTS 


Cfc) 


Fig. l.lj RLC Filter Network and the Corresponding NASAP Input Instructions 



possible set of number assignments for the RLC filter' is shown in Fig. .1.1(b) 

- Overall Appearance of Input Instructions 

Figure 1.2 shows the general appearance of the input instructions. The 
instructions fall into three different categories: a) Control Statements, 
b) Network Element Table, and c) Output Specifications. Each type of ins true 
tions will be explained briefly below. 

Control Statements 

The first line of the set of instruction must start with the word: NASAP 
which simply indicates the beginning of a circuit description. The letter N 
must be in the first column of the card. The word NASAP may be followed by 
any identifying information but only the first 50 columns of the card will 
be printed by the computer. 

The end of the network element table is designated by the ^control card ' 

' " "A 

with the word: OUTPUT -on it. Following the OUTPUT control card are the out- 
put specifications. The final control card appearing after the output speci- 
fications is the card with the word: EXECUTION 1 . 

This card signals the end of the description for the problem and starts 
the computer on the task of analysis. If it is desired to analyze several 
problems in one run, the EXECUTION card can be followed by another NASAP 
card, etc. 

Network Element Table 
General Format 

The network elements and topology are described in a five column network 
element table. (Note that the example in Fig. 1.2 has only four columns 
since no dependent generators are used.) A brief description of each 
column follows ; 
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Column 1: Identifies the element itself., i.e., VI, R3, etc. 

Column 2 and 3: Define the nodes between which the element appears. The 
positive current direction is always from the first to the 
second node. Thus, for passive elements the first node is 
assumed positive with respect to the second, and for active 
elements the second node is assumed positive with respect to 
the first . 

Column 4: Contains the value of the circuit component. 

Column 5: Used only for dependent sources and will he described in 

detail later. 

input Voltage or Current 

The first line after the MSAP control statement contains a description - 
of the applied voltage or current. For example, the statement; VI 1~2~1.0‘ in; 
Fig./oi2 indicates ‘that VI is a voltage of magnitude 1.0 applied between 
nodes 1 and 2 with node 2 being positive. The fourth column contains the 
value of the generated current or voltage and is usually set to 1.0. If 
any other value is used, it will act as a constant multiplier for the transfer 
function under consideration. 

Resistors, Capacitors and Inductors 

Resistors, capacitors and inductors are identified in column 1 by the 
letters R, C, and L, respectively. The letter is immediately followed by 
a one or two digit integer indicating the element number. The fourth column 
contains the component value. The value is a number followed by up to two 
letters which denote a multiplying factor according to the convention given 

in Table fill. The network element table in Fig. li.!2, follows directly from 

v/ _ _ j j 

the circuit diagram of :Fig. £l (b). 
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Table 1.1. Units Following the Component Values 

i e * 


Letters Used 
in NASAP 

Electrical 

Units 

Multiplying 

Factor 

• -Ho letter 

Ohms 

1 

K- 

Kilohms 

103 n 

M' 

Megohms'^ 

los-];- 

F 

Farads 


UF 

Microfarads 

10" 6 

IF 

Picofarads 

10- is 

H 

Henries 

1" 

MR 

Millihenries 

lQ r 3 

UH 

Microhenries 

10-'® 


Output Specifications 

Once the electrical circuit model of the control system has been described, 
we need to specify what we want calculated. This specification always starts 
with a transfer function. 

- Transfer Function 

The transfer functions are always specified as ajratio, ,of- the- voltages?' 
across >or current through a circuit element at the network output over the 
input yoltage or current. The input voltage or current are specified by the 
letters W1 or III respectively and either one must appear in the denomina- 
tor of the transfer function specification. For example, VR3/W1 used in • 

• Fig,..- 1 1^2 implies that we want to calculate the ratio of the voltage across 
R3 to the input voltage VI, i.e., the open circuit transfer voltage function. 
Since R3 was defined by the statement 
R3 4 1 0.8333 

node 4 is assumed to be positive with respect to node 1. If R3 has been 
defined by 
R3 1 4 0.8333 

the sign of the transfer function would be negative . 
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When the program is executed using the circuit description shown in 
Fig. 1.2 and the transfer function specification, the computer generates the 
output shown in Fig. 1.3. The first information printed under the heading 
"NUMBER OF LOOPS PER ORDER" describes the complexity of the circuit -fLowgraph by 
indicating; the number of first order loops, the number of' second order loops, 
etc. This information is only of theoretical interest and can usually be 
ignored during design. However, since the program prints this data, it will 
always be included for completeness. 

Dependent Current and Voltage Sources 

NASAP permits use of any of the four possible types of dependent generators 
(also known as controlled sources): 

a) voltage controlled voltage sources (VCVS ) 

b) current controlled current sources (ICIS) 

c) voltage controlled current sources (VCIS ) 

d) current controlled voltage sources (ICVS) 

Each dependent generator involves two pairs of nodes and two elements: the 
dependent generator itself and a passive element which defines the controlling 
voltage or current. An example of each type of dependent source taken from 
[MO lj is shown in Fig. 1.4. • The node numbers are chosen to indicate that 
they are elements in a larger circuit and the dependent source is indicated 
by a diamond to distinguish it from an independent source which will be 
represented by a circle. 

A description of the NASAP coding for dependent generators is in order. 

The VCVS shown in Fig. 1.4(a) is specified in NASAP by the instructions 
found immediately below the figure. The first instruction specifies a de- 
pendent voltage source, V3, from node 8 to node 3 (node 3 positive). The 
value of the source is 867.3 times the voltage across capacitor C2. The 



NUMBER Of loops Per order 


1= 5 

2 = 4 


TRANSFER FUNCTION VR3/VV1 


( 

( 5 « OOE 00 

HCS)= 3 # 333E-01*--- -- 

C 

( 3 « OOE 00 


+ 1 ♦ OOF 00 S ) 


2 ) 

+3,oaE 00 s + 1 , OOE 00 s ) 


ZEROS of TRaNSFfR function 
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1 “.15000E oi 

2 ** « 1 5000E 01 


FUNCTION 
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Figure 
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V3 8 3 876. 3 VC2 
C2 5 1 2. O'UF 

(a) VCVS 
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II 6 1 50 IR2 
R2 3 1 350 . 

(b) ICIS 



L3 2 4 5. OMH R7 \z 7 10K 

12 6 1 0. 8 VL3 . YE 17 8 100. IR7 

(c) VCIS (d ) ICVS 


Figure a; k, \ 


Examples of the Four Types of 
Dependent Generators and the 
Corresponding NASAP Instructions 



second instruction specifies a capacitor from node 5 . to . node - X, 6f 2;0pf . Thus 
node 5 of the capacitor is assumed positive as noted in Section: Network Element 
Table, page 1.4. The polarity of the passive controlling element will always be 
shown on the circuit diagram to help insure that the nodes will be numbered in 
the correct sequence since the passive element may appear far from the controlle 
source in the network element table. 

Similarly the ICIS shown in Fig. 1.4(b) is specified by the instructions 
below it. That is, a current source, II, from node 6 to node 1 with a value of 
50 times the current flowing through R2. R2 is defined so that node 3 is assume 
positive and thus the current is assumed to flow from node 3 to node 1. 

The corresponding NASAP instructions for the VC IS shown in Fig. 1.4(c) 
and the ICVS shown in Fig. 1.4(d) are written to indicate that it does not 
matter whether the dependent generator or the passive controlling element is 
specified first in the network element table. 

Sensitivity 

The MSAP program is capable of determining the sensitivity of the transfei 
function to changes in the value of any single element. The request for a sen- 
sitivity analysis is optional and, if it is desired, is specified following the 
transfer function specification. This has the typical form: 

VR11/W1/L3 

This output specification requests the sensitivity of the transfer function 
VRll/wi with respect to element L3. 

The frequency range over which the sensitivity is analyzed is the same as 
that for which the frequency response is calculated so they are both defined 
in the same output specification. For the RLC filter of Fig. 1.1, Moe [MO l] 
tabulates and displays the response for FREQ -1*2 I'. O' 0.0'5-. 



Further Discussion of Dependent Sources 


As stated in Gaertner's Coding manual, dependent sources can be related 
only to the voltages or currents of passive elements of the flowgraph. How- 
ever, the controlling current or voltage in the circuit may be that of an 
active element. Such a situation is illustrated in Fig,.- 1*5 ' by the circuit 
realization for a normalized three-pole low-pass Chebyshev transfer function. 
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Desired Transfer Function: VC3/W1 

Fig. 1.5 

This circuit is an application of a general circuit introduced by Cooper 
and Harbour t [CO-l] . ~ Note that the controlling voltage, of the dependent 
source V2, is the. independent voltage source VI. However if the following 
card 

V2 1 6 -2 W1 

where = - 2 
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is included in the input list, the error message 
INPUT CODING ERROR IN COLUMN 12 

will result. The error is the second letter V since the program in only 
searching for the letter R, L, or C. 

However, this difficulty is easily resolved by connecting a resistor R 6 
across the terminals of VI (i.e., between nodes 2 and l) and by making the voltage 
source' dependent on the current through ^R 6. The current through R 6 (IR6) equals 
V1/R6. 

Thus VI - (R6)IR6. Since V2 - 1^ VI, then 
V 2 = \ * IR6 

Thus the dependency value between V2 and JR 6 is .'r 6. Note that R 6 can have 
any nonzero numerical value . If, as is the case in this example, R 6 is chosen 
to have a value of 4 ohms, input list must include the following two cards; 


R6 2 1 4 
V2 1 6 -8 ir6 


(see input list in Fig. 1.8) 


The presence of this extra resistor heither increases the complexity of the 
flowgraph by generating additional loop sets nor does it affect any of the elec- 
trical properties of the original network. 

This can be seen by comparison of -Figs*. .3^6"' ' and - -Fi'gv 1 . 6 is* a partial 
flowgraph for the original network in Fig. X.5T • ;Eig.. i -7 'is a partial ;flow- 
graph for the addition of r 6 to the network. 

Note that in both Figa.. 176 ^ahd” l>7~f^ere is only one path from the voltage 
node of VI to the voltage node of V2, 

The computer results for the circuit in Fig.,- 1 . 5 r are -Jiveh' in. Figure" 1,8."' 

f~ - - ' 4 . ‘ ! . 

A similar procedure is followed when a dependent source is a function of the 

current of a current source. In this case a resistor R is added in series with., ' 

s - * 





the current source. The dependent source is then made a function of the voltage 

across this resistor — the dependency value being k /R where k is the dependency 

I s i. 

value between the dependent source and the current source. As in the previous 
case addition of this resistor to the network under consideration, will in no 
way affect the flow graph or the electrical properties of the original network. 
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iB ERROR messMs ' 

- - * - — - — . vi*?:. - - 

Here are provided a few comments on the meaning of selected error messages 
available in NASAP. 

In subroutine GETCON, there are two error messages. 

The message 

ALGORITHM FAILURE REVIEW CIRCUIT CODING RULES 
usually occurs when the -NASAP program is unable to select a legitimate tree from 
the available resistors or capacitors. This case exists when the candidate type 
1 or 2 elements form an R or C loop or when all elements connected to a node are 
type 6 or 7 elements; i.e., these excluded elements from a*;lcut - set . If r for' 
example, the transfer function VR4/W1 is desired for the circuit in ;Fig . 1;9« 


C2 C3 



then the above error message will be printed out. This is necessitated by the loop 
formed by the elements V2 and r 4' which must be in the tree . This difficulty can 
easily be resolved by calling for the transfer function W2/W1 which is identical 
to VR4/W1. 

The other error message is 

FLOW GRAPH HAS TOO MANY CONNECTIONS EXECUTION TERMINATED 


JL7— 



The connections of the primitive flowgraph are stored in the two dimensional con- 
nection matrix LCONC. The size of this matrix is such that it is assumed that 
no more than nine (edges) branches will emanate from a single flowgraph node. Thus 
if ten or more branches must emanate from a node, the above error message will be 
printed and execution, stopped. This limitation is completely arbitrary and exists 
to restrict the size of the matrix LCONC. If one wishes to increase the maximum 
number of node connections, then the number 10 which appears twice on lines ^010, 
5460, and 6220 and which appears once on line 5110 should be changed to a suitable 
integer N which -is defined as'the maximum number of node connections (m) plus one. 
Also the numeral 9 on line 6250 should be changed to the integer M. When these 
changes are made, the above error message will result when the number of node con- 
nections exceeds M instead of when the number of connections exceeds nine. 

In subroutine FSORL, there is the error message 
FLOWGRAPH FIRST ORDER LOOPS EXCEED 927. 

This limitation on the actual loops of the flowgraph, the so-called first order 
loops, is completely arbitrary. It is necessary to store the first order loops 
since the higher order loops are determined from these first order loops. However, 
if this error message is printed out, it usually means that the tree selected was. 
a rather poor choice. A more complete explanation of tree selection to minimize the 
number of loops of a flowgraph will be given later. 

In subroutine BODE the error message is 

PROGRAM RESTRICTS GRAPH TO 250 STEPS . 

The sub script for the evaluation of the transfer function at various frequencies 
cannot exceed 250. Thus, if the information on the FREQ input card results in 
more than 250 frequency evaluations, the NASAP program will not perform any of the 
calculations included in subroutine BODE including the sensitivity calculations 
(if this has been called for). The program will jump to the time response calcula- 
tion if this has been called for. Otherwise, .execution will cease after the above 


error message. 
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In subroutine SENS there is the error message 

SENSITIVITY PLOTS BEffinSCTED TO 120 POINTS . 

If the input information on the FREQ card results in more than 120 frequency 
points , the sensitivity results in tabular form will be printed out. However the 
3-curve plot and the phase sensitivity plot will be deleted. 

At present, extensive diagnostics and debugging capabilities are lac Icing . Input 
data can be written on a field- free format, and -a circuit tree selection ; need not 

be specified by the user since it is done .internally^ by the program. As is some- 

~~~ "W 

times the case, when a proper, tree cannot be found, a message is 'printed out indi- 
cating the difficulty ‘and * ■(there it. occuFs71.n jthe/circuit . . Detailed 7 discussion- ofl 

'“t — * * ajHMf ' 

the tree selection algorithm is found in Chapter II. 
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IC DOCUMENTATION OFJNASAP 69/1 VERSION USED AT THE UNIVERSITY OF PENNSYLVANIA 
Summary 

This section details the differences in input data and output results of the 
University of Pennsylvania version as compared to the standard NASAP 69/1 as des- 
cribed in the Gaertner coding manual. [GA-l] 


Input Cards 


The NASAP version used at the Moore School of Electrical Engineering in the 
University of Pennsylvania (hereafter called the MSE-NASAP version) requires the 
user to supply two additional data cards immediately after the first data card 
(the NASAP PROBLEM card) . The first of these cards indicates the frequency units 
to be used in the evaluation of the transfer function. There are four permissible 
entries for this card (starting in column l) 


RADIANS 

HERTZ 

CYCLES PER SECOND 
NONE 

Only the first two letters (those underlined) must be correct since the NASAP 
program checks only these letters. Any other information can be included on this 1 card 
since the program evaluates only the information contained in ’columns 1 and 2. 

The second card indicates the type of time response desired by the user. There 
are four choices (starting in column l) 

IMPULSE RESPONSE 
STEP RESPONSE 
RAMP RESPONSE 
NONE 

The programs again only evaluate the data in columns 1 and 2. 


20 



Table 1.-2' shows the input listing necessary on the MSE-NASAP version to find 
the step response for the circuit of Fig. 3Ti.107'' The response is the voltage ^across 
Cl and the excitation is a step voltage with a magnitude of 3 volts. The BODE 
plots of the specified transfer function are not being requested. 


VI 


r—* VW 

o 

R1 




> ^ 

; 0.5F 

Cl 

0 


Transfer Function: V C1 /W1 
Fig. 1.10 


Table' 1.-2. 

NASAP PROBLEM 
NOME 

STEP RESPONSE 
VI 1 2 3. 

HI 2 3 1 
Cl 3 1 0.5F 
OUTPUT 
VCl/Wl 
TIME 2.0 
EXECUTE 
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Printed Output Results 


The MSE-ISASAP version prints out much data on the primitive flowgraph which 
is used within the computer to determine the user- specified transfer function. 

The additional output data is printed between the printed listing of the input 
data cards and the NUMBER OF LOOPS PER ORDER table— both of which are printed by 
the standard NASAP package. 

As an example. Fig. 1.11 gives printed output of a legitimate input listing 
for circuit 6 introduced in Hutton* s report. [HU 1] Fig. 1.12 gives the 5 com- 
puter sheets of additional data printed out by the MSE-HASAP version based on . 
the listing of (Fig. /l. 11. : 0n the sheet immediately^ following^tHe^input list - sheet--, 
are printed two matrices. The first matrix is the compressed cut set matrix. 

This matrix has (b + 2) rows and (& + 2) columns where 

b ss number of branches in the tree = n - 1 
n = number of nodes in circuit 
- & = number of links in the . co-tree 

It must be recalled that 

b + X — X ( 

where X is the number of elements in the circuit (For the circuit described in 
Fig. 1.11, X - 30 and n = 15). This matrix is useful in that it shows which 
circuit elements have been selected as tree branches by the NASAP tree- selection 
algorithm. 

As each circuit element is inputted, it is assigned an integer flag (be- 
ginning with ('unity J>. Thus for the circuit listing of Fig. 1.11. 
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Element 


Integer Flag 


VI 1 

r;6 “ 2 

E15 3 

V3 30 

An element in the compr e s se m ouqse u , matrix can have one of three values -1, 

0, or +1. The top and "bottom rows and the left-most and right-most columns of 
the compressed^ cut set i matrix are used for identification purposes described in 
the following paragraphs. 

The non-zero entries in the left-most column are the integer flags of those 
circuit elements which have been selected as tree branches. The corresponding 
entry in the right -most column gives the tree hierarchy value of the circuit 
element (see description of tree -selection algorithm). Hence in Fig - -we see : 

that VI (flag s l) and R 7 (flag - 9) are tree branches and both have a hierarchy 
value of 2. 

Similarly the non- zero entries in the top row are the integer flags of the 
circuit elements that are links of the co-tree. The bottom row gives the hierarchy 
value of these elements. For instance, HI (flag = if) and R2(flag = 5) are links 
of the co-tree and have hierarchy values of if. 

The information contained in the compre s seel cut se t^matr ix is used to develop the 
flowgraph connection matrix which is a mathematical description of the primitive 
flowgraph. This matrix consists of X (= number of elements) rows. The entries in ~ 
the left-most column represent either the voltage node or the current node of 
the circuit elements (depending upon whether the elements are branches or links 
respectively) identified by the integer flags numbered consecutively from 1 to X. 
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If the left-most entry of a particular row of the flowgraph connection matrix 
represents a voltage node of an element, then each of the other entries of this row 
indicates a connection from the voltage node of the left-most, entry to the voltage 
node of the corresponding entry. There is one exception— if one of these entries 
refers to a voltage controlled current source, then there is a connection from the 
voltage node of the leftmost entry to the current node of this particular entry. 
Conversely, if the left-most entry of a particular row of the flowgraph connection 
matrix represents a current node of an element, then each of the other entries of 
this row indicates a connection from the current node of the left-most entry to 
the current node of the corresponding entry. There is one exception — if one of 
these entries refers to a current controlled voltage source, then there is a con- 
nection from the current node of the left-most entry to the voltage node of this 
particular entry. A few examples from Figure.. 1.12 will ciari^-thisf.dfescf iptio^n^ iThe 
second row of the flowgraph connection matrix in Figure ;£l2 contains 

2 -8 -11 0 

From the cutset matrix it is seen that the element with integer flag 2 (i.e. R6) is 
a tree branch. Thus there are two connections (with a value -1) from the voltage > 
node of r 6 - one to the voltage node of R5(flag = 8) and the other to the voltage 
node of R9( : ?lag =11). The zero indicates the end of the row. 

As a second example, the fifth row of this connection matrix contains 

5 9 -19 20 0 

Since R2(flag = 5) is a link in the co-tree, there is a connection with a value 
of +1 from current node of R2 to the current node of R7(flag = 9). Also there is 
a connection (value = -l) to the current node of Cl (flag = 19) and one (with 
value = -1) to the current node of C2(flag = 20). Use of Ohm’s Law yields the 
passive element (R,L,C) joining the voltage and current nodes, since the voltage 
across an element is related to the current flowing through the element by means 
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of the element's impedance. 

The next printed output consists of all the closed loops (first order loops) 
contained in the flowgraph determined "by means of the flowgraph connection matrix. 
The heading of this output block of data is 

FIRST ORDER LOOPS BY CONSECUTIVE LOOPS 

Each loop is defined by the nodes contained in the loop. This output shows the 
order in which the nodes have been found by NASAP for each loop. Each loop is 
given an identification integer starting with unity. This identification integer 
is shown in the left-most columns. The remaining integers in a given row refer 

to the integer flags of the circuit elements. 

■ “•* ’ * • 

In Figure. 1 .^2jwe see tha t: the -flowgraph. .under discussion^ has. ,a total of 102' 
first order loops. Let us examine one of these loops more closely — say loop 
81 defined by 

-6 2 6 7 -6 

The first and last integers are always identical since it is the starting point 
of the loop. From the input listing we have 

Element Integer Flag 

R3 6 

II 2 6 

Bk 7 

Each loop is found by the path-finding procedure on the flowgraph connection 

matrix (for details see the Pot ash-Mc Name e User's Manual). From the .cutset 
matrix we note that R4 is a tree branch while II and R3 are co-tree links. 

To illustrate this we now show how this loop 8l can be found from the flow- ' 
graph connection matrix. Start at the current node of R3 (R3 is a link). There 
is a connection from this node to the current node of the current source II (i.e., 
an integer 2 6 in row 6 of the connection matrix).. There is a connection from 
the current node of II to the current node of r 4 (i.e., an integer 7 in the 26th 



row of the connection matrix) . Since r4 is a passive element, there is a con- 
nection (with a value = k4) from the current node of r4 to the voltage node of 
R4 (recall that r 4 is a tree branch). From the 7th row of the connection matrix, 
(integer flag of r4 = if) we observe that £6 Vis an entry. Thus there is a connec- 
tion from the voltage node of R4 to the voltage node of R3 (with a value of -l).. 
Since R3 is a passive element that is a co-tree link, there is a connection 
(with value = ^3) from the voltage node of R3 to the current node of R3 (the 
starting point)— thus completing the loop. 

The higher order loops, the sets of non- touching first order loops, are 
easily found by assigning each first order loop an integer value based upon the 
nodes contained in the loop (each node is identified by the circuit element 
integer flag). This integer value is stored in the one -dimensional array LOOP 
defined by 

|N. | , 

LOOP(J) = t 2 * ± 

k=l 


where J is the loop identification integer 

refers to the kth flowgraph node in the Jth loop 
n = number of flowgraph nodes in the Jth loop 
As an example, for the loop numbered 79 in FigT'£^l2, we have j 


|-6|-i 1-7 l.'-l 

LOOP (79) =2 + 2 r 


2 5 + ^ 


loop(79) 


(96), 


'10 

= (110000) 2 

(Note: see Potash-McNamee Manual for details of how higher-order loop are found 
by use'of_the array LOOP). 

Finally the sheet with the heading 


ORDER OF LOOP 


LOOP NUMBERS 



is printed out to assifcythe user in locating the point of termination if the 
number of flowgraph loops is so large that -the allowable computer time is used 
up before completion of the analysis. 

The information of this sheet is presented in this manner: 

1. Every 50th second order loop is printed out with the identification 
integers of the first order loop comprising the second order loop. 

2. When the number of some jth order loop (j = 2,3,4,...) equals 500, 
the number of all loops of all order at this point of the loop 
enumeration procedure. (Note: see the Potash-McNamee Manual for 
details of the HIGORL subroutine which determines the higher order 
loops from the first order loops ) . 

From Figure^ 1.12 we see that the output data 

2 18 88 

is the fourth of this type in the output data. Consequently the second order 
loop consisting of loop 18 and loop 88 is the 200th second order loop found by 
the NASAP subroutine called HIGORL. (Note: from the list of first order loops, 
loop 88 consists of nodes 10 and 12 which do not appear as nodes in loop l8) . 
Immediately following this line of output, we observe that the number of 9th 
order loops found by NASAP equals 500. When this occurred, NASAP had found 209 
2nd order loops, 987 3rd order loops, etc. 

It should be noted that, if a flowgraph has less than 50 second order loops 
and the number of each of the nth (n = 3,4,5, .*.) order loops of this flow- 
graph is less than 500, only the heading 

ORDER OF LOOP LOOP NUMBERS 

will be printed in the MSE-NASAP version. 

In order to find all the higher order loops and at the same time to avoid 
repeating any of these loops, the HIGORL subroutine selects each first order 
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loop in order of its identification integer and determines all the higher order 
loops formed by this selected first order loop and those loops having identification 
integers greater than that of the selected loop (see Potash-McNamee manual for 
details ) . 

With this procedure in mind, we can use the higher-order loop data printed out 
by the MSE-NASAP version to determine approximately where the HIGORL subroutine 
was in the higher order loop finding process if execution is terminated before the 
desired transfer function is obtained. This will enable the user to determine 
whether the problem should be executed with a longer running time specified or 
should be cancelled since the generated flowgraph is too complicated (with 
regard to running time) for the computer in use. 

As an example, let us assume that the printed output of Figure ; 1.12. endedwith 
the line 

2 68 88 

There was no sheet with the heading 

NUMBER OF LOOPS PER ORDER 

and no sheet listing the specified transfer function. From the list of first 
order loops, we see that this flowgraph contains 102 loops. From the above 
output line, we see that the loop mnder selection by HIGORL is the loop num- 
bered 68. In other words, when this line was printed loops 69 through 101 
have not yet been selected as the starting loop in the higher order loop 
finding process. 

Since the number of higher order loops generated by the starting loop de- 
creases as the identification integer of the starting loop increases, knowledge 
that loop 68 (in a group of 102 loops) is the present starting loop enables the* 
user to decide whether or not to rerun the problem. It must be noted that the 
last line printed does not mean that this was the last loop set found by the 
HIGORL subroutine before termination . It is possible that as many as b9 more 
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2nd order loops and many more nth (n = 3,4, 5 * . . ) order loops have been found 
since the printing of the last output line and before the termination. 

Following this output data, the output results of the MSE.-HASAP version con- 
form with the output data of NASAP 69 /I. One remaining minor difference is the 
order of the sensitivity output data. 

The order for NASAP 69/1 is: 

Table of Sensitivities of ReH, ImH, |h|, Phase H as a function 
of frequency. 

2. Table of the Logarithm ^of the above sensitivities as a function 
of frequency, 

3. Plot of 3 Sensitivity Expressions. 

4. Plot of the Logarithm of the Sensitivity. 

5* Table of the Sensitivity Function as a function of frequency. 

6 . Plot of the Logarithm of the Absolute Value of the Sensitivity 
function. 

7. Plot of the Phase of the Sensitivity Function. 

8. Pole and Zero Sensitivities. 

On the. other hand, the order of the MSE-MSAP version is; 

1. Table of Sensitivity Function as a function of frequency. 

2. Plot of Logarithm of Absolute Value of Sensitivity Function. 

3. Plot of Phase of Sensitivity Function. 

4. Table of Sensitivities of ReH, ImH, |h|. Phase H as a function 
of frequency. 

5. Table of Logarithms of the above sensitivities as a function 
of frequency 

6 . Plot of 3 Sensitivity Expressions. 

7. Plot of the Logarithm of the Sensitivity of the Phase H. 

8. Pole and Zero Sensitivities. 
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CHAPTER IX 


NASAP TREE SELECTION ALGORITHM— USER OPTIONS 
IIA : C^NERAh' DESCRIPTION ' 

Although the graph representation of an electric network usually has a large 
number of possible trees (i.e., a structure .containing n-1 branches which inter- 
connect the n nodes of the circuit without forming any closed paths), NASAP has 
an algorithm that selects only a particular tree configuration. This tree is 
the basis for subsequent circuit analysis. 

Each electrical circuit element is assigned a Type number as shown. 


Independent voltage source 1 
Dependent voltage source 2 
Capacitor 3 
Resistor 4 
Inductor 5 
Dependent current source 6 
Independent current source 7 


Elements of type 1 and 2 are always included in the tree while type 6 and 7 ele- 
ments are never included in branches of a tree . If there are not enough elements 
of type 1 and 2 to form a tree, a search is made of type 3 elements (i.e. capa- 
citors) starting with the first capacitance listed in the input data and working 
down the input list. If a tree is .still not found after searching through all the 
type 3 elements, a similar search is made of all type 4 elements (selecting 
those type 4 elements that do not form closed paths and neglecting those that 
do). If a tree does not result, a search is made of type 5 elements. If a 
tree is not found after this; search, an error message will be printed. 

The element type categories 2 and 6 need a further explanation. The elements 
in category 2 include not only dependent voltage sources, but also those elements 
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whose voltages control the voltage or current of some dependent source or .whose 

f , 

voltage is the required output variable or * input variable. For example, the 
following is a legal NASAP input record. 

13 3 4 5*2 VR3 

Since the voltage across resistor R3 controls the dependent current source 13, 
resistor R3 will be assigned element type 2 not element type 4. Also the fol- 
lowing is a legal NASAP output record 

VL5/W1/R2 

Since the voltage across inductor L5 is the desired output variable. Element 
type 2 not element type 5 will be assigned to inductor L5 . Similarly since VI is 
the input, it will be assigned element type 2 not element type 1. 

Similarly the elements in category 6 include not only dependent current sources 
but also those elements whose currents control the voltage or, current of some de- 
pendent source or whose current is the required output variable or input variable. 
In 

13 3 4 5.2 XR3 

Resistor R3 is assigned element type 6 since the current through R3 controls the 
dependent current source 13. In 

IC5/II2 

napaaitor C5 .is assigned element type 6 since the current through C5 is the re- 
quired output current. Also, 12 will be a type { 6 -' element (not type 7) since it 
is the input variable. 

Due to the search ^technique 7 that" reads down ^he : input ^iis'tTl^Hh^for^lementsi. 
to be branches of a tree, the actual tree selected by NASAP can be varied simply 
by rearranging the order in which the elements are placed in the input list (NOTE: 
There is no restriction as to the order in which the elements are listed in 
NASAP). However, not all of the possible tree configurations can be selected 
by NASAP due to the requirement that all type 1 and 2 elements must be included in 
the tree while all type 6 and 7 elements cannot be included in the tree . 

37 



I IB ILLUSTRATIVE EXAMPLE 


The complexity of the flowgrapb is greatly dependent on the particular tree . 
By complexity Is meant the number of loop sets present in the flowgrapb. The 
values of these loop sets are used In HASAP to calculate the transfer function by 
means of Mason’s formula as extended by Rapp for the closed signal -flow graph* 
Since much time is consumed by the NASA? algorithms in finding the loop sets, 
selection of the tree that minimizes the number of loop set can save computer 
time as well as increase the accuracy of the coefficients of the polynomials in 
the transfer function. 

In Fig. 2-1 i 0 shown a transistor with known h -parameters * The input resis- 
tance of the circuit is to be calculated. 



Fig. 2.1 4 A Common Emitter Transistor 


The RASAP equivalent circuit model is given in Fig, 2*2, 


® B 


R2 © 


® C 


XI® 


-1M IM 



R I 



T *E 

!k£R4 


R3 

40k 





Fig . 2.2. Equivalent Circuit for Transistor in Fig. 2-1. 
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II is the independent current source. R1 and r 6 are positive and negative re- 
sistances of equal numberical value. These resistances, when added to the ■ 
network, yield an element whose voltage equals the input voltage and at the same 
time does not load the circuit (the parallel connection of a positive and nega- 
tive resistance yields an infinite resistance). Thus the required transfer 
function for the input impedance is 

VR1/IH 

The network in Rig. 2.2 has five nodes; thus the tree for this circuit must have 
four elements . The elements Rl, VI, and R3 are type 2 elements and must be in- 
cluded in the tree while II, 12, and R2 are type 6 elements and cannot be part 
of the tree. Thus one more element must be selected to form the tree. The 
partially completed tree structure is shown in Fig. 2.3. 



r 6 cannot be used as a tree element since it forms a closed path with Rl which is 
an element in tree . However either R5 or R 4 can be used to form a tree . Thus the 
HASAP tree selection algorithm will pick that resistance which is listed first in 
the input data. Thus if the input data is listed as follows, 

n i 2 i. 

12 5 4 100 IR2 

VI 4 3 0.0001 VR3 
Rl 2 1 1M 
R2 2 3 IK 
R3 5 4 40K 
R5 5 1 10k 
r 4 4 l ik 

R6 2 1 -1M 

OUTPUT 

vRi/m 














NA'SAP I IT PROBLEM 

ohXL ~ M /ts^ivU<3l zs> 


NONE ' 
NONE 


11 1 2 
1-2 5 A 

1. 

100 IR2 

Vi 4 3 

0.0001 VR3 

R1 2 1 

XM 

R2 2 3 

IK 

R3 5 4 

40K 

R5 5 1 

10K 

R4 4 1 

IK 

R6 2 1 
OUTPUT 

-XM 

VR1/II1 


EXECUTE 



NUMBER OF LOOPS PER ORDER 


1= i4 
2= 29 

3 - 10 


TRANSFER FUNCTION VR3/II1 


( ) 
( l.OGE 00 ) 


H ( S ) " 8.03-3E 0*-* 


I l 

{ 1 • OOE 00 ) 


ZERO OF TRANSFER FUNCTION 


NONE 


PO.LE OF TRANSFER FUNCTION 
NONE 


Fig, g.j? 
hi 




R5 "will be included in the tree elements. The chosen tree is shown in Rig. 2.4. 


3 <D © © 


Ri \^-^R5 

0 _ _ Rig. 2.4 

'The output results are given in Rig^T^n?^ - „,. r 
However, if the input listing format is given as 


11 

1 

2 

l. 


12 

5 

4 

100 IR2 


VI 

4 

3 

0.0001 VR3 


R1 

2 

l 

IM 


R2 

2 

3 

IK 


R3 

5 

4 

40k 


r4 

4 

1 

IK 

(Note that R5 is now 

r 6 

2 

i 

-1M 

the last element des- 

R5 

5 

l 

10K 

cribed) 

OUTPUT 




TOi/m 




/■cne^tree selected by Ml SAP will now contain R4 as shown in Rig,~'~2.6. 


0 © © 0 



The output results are shown in Rig. 2,7* 

Note that both coding formats yield the same transfer function, as expected. 
However, the tree with R4 yields a flowgraph with 29 loop sets as compared with 
53 loop sets in the flowgraph obtained from the tree with R5. Thus the number 
of loop sets is almost halve d^by the proper choice of a tree for this circuit 
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configuration . This reduction of loop sets is apparent in the execution time. 
The tree with R5, KASAP required an execution time of 24.3 seconds on the RCA 
Spectra 70/46 while the execution time for the tree with r 4 was 22.6 seconds, a 
saving of 1.7 seconds. 


43 



IIC THE OPTIMUM TREE 


There is a procedure to select the optimum tree (optimum in the sense that 
such a tree will minimize the total number of loops in the primitive flowgraph) 
that does not require a great deal of effort on the part of the NASAP user. 

The details and proof of this procedure are given in reference [Z0-l3. The pro- 
cedure begins by putting all type 1 and 2 elements in the tree. Next all capa- 
citors (type 3 elements) are included in the tree. If a capacitor forms a loop 
with some type 1 or 2 elements, it is removed from the tree. If two capacitors 
form a loop with some other tree branches, then the user arbitrarily picks one of 
these capacitors to be a tree branch — realizing that the capacitor picked to be 
the tree branch must precede the other capacitor in the input lists. All type 
6 and 7 elements and any element that forms a loop with the chosen tree branches 
are put in the co-tree. Note that if the sum of the type 1, 2, and 3 elements 
equals the number of nodes minus one and if these elements do not form any loops, 
then the NASAP tree selection algorithm will have picked a tree after a search 
of all capacitors. The user thus will not be able to vary the tree. 

If, on the other hand, the number of resistors in the circuit is greater than 
the number of elements necessary to complete the tree, the NASAP user will have 
some flexibility in the type of tree selected by NASAP by permitting the resistor 
input cards. 

Once the user has selected a tree, the optimum tree search procedure goes 
as follows: 

Each link forms a loop with some of the branches of the tree. For each link, 
the number of tree branches in each loop is recorded as well as the specific 
branches that form the loop. The Branch Count is the sum of the number of tree 
branches in each loop. The Circuit Count is the sum of the number of loops for 
a specified tree branch. The Branch Count will always equal the Circuit Count. 
However the tree that yields a smaller Branch Count will yield a flowgraph with 
fewer loops. 
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As an example let us reconsider the input impedance circuit given above. 


Recall that it was shown that VI, HI, and R3 must he included in the tree while 
II, 12, R2 and R 6 must he links in the co-tree. However either r 4 or R5 can he 
in the tree depending on the order of the input cards. If r 4 is a branch of the 
tree, then there are five links in the co-tree - H, 12, R2, R 6, and R5. 

The link II forms a loop with branch R1 
The link 12 forms a loop with branch R3 
The link R2 forms a loop with branches Rl, VI, R4 
The link R 6 forms a loop with branch Rl 
Finally, the link R5 forms a loop with branches R3, R 4. Table 2.1 .summarizes 
this information. ("Branch" refers to a tree branch.) 


Links 

Branches/Loop 

Branches 

Loops/Branch 


Table 2.1 

II 12 R2 

1 ' 1 . ,3 

VI Rl 

l 3 


R6 R5 '^branch count 
* / _2 . 8 
R3 R4 * circuit count 

2 2 8 


Note that the Branch Count equals the Circuit Count, as required. 
However, if R5 is made a tree branch, then the co-tree links are ZL, 12, 
R2, r6, r4. 


The 

The 

The 

The 

The 

Table 2.2 summarizes this 


link II forms 
link 12 forms 
link R2 forms 
link R6 forms 
link R4 forms 
information. 


a loop with branch Rl 
a loop with branch R3 
a loop with branches Rl, VI, R3, and R5 
a loop with branch Rl 
a loop with branches R3 and R5 
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Table 2*2 


Links 


11 12 B2 ft 6 RU Branch count 


Branche s/Loop 


1 1 1 U 


2 9 


Branches 


VI R1 R3 R5 Circuit count 



1 3 3 a 9 


Both the Branch Count and Circuit Count equal nine. Since the Branch Count 
for the tree containing r 4 is less than that for the tree with R5, then the number 
of loops for the flowgraph formed from the tree with r 4 will be less than that 
obtained from the tree with R5* As noted above, there are 29 loop sets including 
9 first order loops in the flowgraph formed from the tree as opposed to the 53 loop 
sets including lh first order loops in the flowgraph formed from the R5 tree. 

Illustrative Example 

The circuit shown in Fig, 2*8 from iMO-i] illustrates the case when different 
trees yield the same branch count. 



X I * 500 V R £ 

Desired transfer function VR0/W1 


Fig. 2*8* An EOG Filter from [MO-l] 


It is apparent that VI, R6, and r 8 (type 2 elements) and Cl must be branches 
of the tree while II, R5, and R7 must be links in the co-tree since they form loop 
with the above mentioned tree branches. For the "six node circuit in Fig. 2.8 a 
tree has five branches. Since four of these branches have already been designated 
the last branch can be chosen from resistor Rl, R2, R3 or r 4. Thus four different 
trees can be formed by the HASAP algorithm. The question remains which of these 
four trees will yield the fewest loops in the flowgraph. 

If R2 is selected as the fifth tree branch, then the links of the co-tree 
are II, R5, R7, Rl, R3, and r 4. Table. 273; gives the pertinent information on the 
number of branches in 'each loop formed by each link and the number of loops that 
each branch is a part of. 


Table 2.3: Tree with R2 


Links 

11 

R5 

R7 

Rl 

R3 

Branches/Loop 

1 

3 

1 

3 

4 

Branches 

VI 

R 6 

R8 

Cl 

R2 

Loops /Branch^; - 

2 

4 

4 

4 

- 3 


The Branch Count for the R2 tree is 17. 

■’* If "instead we let r 4 be a branch of the tree; then the links become II, 

R5, R7, Rl, R2, and R3« Table 2.4 gives the corresponding branch and loop in- 
formation. 


Table 2.4: Tree with R4 

Links 

Branches/Loop 
Branches 
Loops/Branch 

Here the Branch Count is 15. Since the branch count for the r 4 tree is 
smaller, then this tree is a better choice than the R2 tree in terms of fewer 
loops for the flowgraph. 


II 

1 

VI 

2 


B5 

3 

R 6 
2 


R7 

1 

r8 

4 


Rl 

3 
Cl 

4 


R2 

5 

r4 

3 


R3 

2 


If we next let R3 be a tree branch, the co-tree links become XI, R5, R7, 

Rl, R2, and R4. Table 2,5 gives the branch and loop data for this case. 

Table 2.5: Tree with R3 

Link II R5 R7 Rl R 2 R 4 

Branches/Loop 13 14 4 2 

Branch VI R 6 r 8 Cl R3 

Loop/Branch 2 2 4 4-3 

Note that the branch count is 15, the same as that obtained for the r 4 
tree. 

Finally Table 2.6 gives the branch and loop data for the case when the tree 
contains the element Rl. 

Table 2.6: Tree with Rl 

Link II R 5 r7 R2 R3 r4 

Branches/Link 13 13 4 3 

Branch VI R 6 r 8 Cl Rl 

Loops/Branch 2 2 4 4 3 

The Branch Count for the Rl tree is also 15. 

Thus three of the four trees have the same low Branch count. Thus the 
other criteria^ the Branch product and the Loop Product, must be used to deter- 
mine which of these three trees will yield the flowgraph with the fewest loops. 

The Loop Product, defined as the product of the number of loops involving each 
branch, is the same for the three trees; namely 192. However, the Branch Product, 
defined as the product of the number of tree branches in each loop formed by each 
co-tree link, is different in each case. The Branch Product is Ix3xlx3x5x2 = 90 
for the R4 tree while the Branch Product for the R3 tree is Ix3xlx4x4x2 = 96 . 

The Branch Product for the Rl tree is Ix3xlx3x4x3 = 108. It is found that the 
tree having the smallest Branch Product (r 4) yields the flowgraph with the fewest 
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Note that the Branch Counts and Branch Products of hoth trees are equal. Since 
the Loop Product of the C2-C5 tree is smaller, it would seem that this tree would 
yield the flowgraph with the fewer loops. 

However, the modified Branch count of the C3-C6 tree is smaller than that of 
the C2-C5 tree. This criterion indicates that the C3-C6 tree gives the fewer flow- 
graph loops. This, in fact, is the case as the computer results of Figs. 2.l8, 

2.19, 2.20, 2.21 show. The tree with C3 and C 6 generates a flowgraph of 969 
loops including 19 first order loops while there are 15 29 loops including 25 
first order loops in the flowgraph formed fromVthe C2-C5 tree. 

From Dunk ^ and Chan,- [DU 1] it has^been shown that the star tree - (aT treeT in-which 
all the branches have a common node) yields the minimum number of flowgraph 
loops. The more star-like the tree structure, the fewer the number of loops in 
the corresponding flowgraph. 

In the circuit of , Pig. 2.17 > all of the type 2 elements are connected together 
at node 1. Thus 5 of the 9 branches are joined at a single node. Examination of 
this circuit reveals that two other elements! (r 4 and r 8) are also connected to 
node 1. If these resistors were branches of the tree, then 7 of the 9 tree 
branches would be connected to a single node — a tree structure that is definitely 
more star-like than the two trees described above. However, since r 4 and R8 
are resistors and type 4 elements, they are not considered for eligibility as 
tree branches by the NASAP tree selection algorithm until all type 3 elements 
(i.e. capacitors) are considered. As is shown above, the capacitors are so 
connected in the circuit that they do form legitimate trees (in fact, 4 trees 
depending upon the input listing) . Thus r 4 and r 8, as type 4 elements, can never 
be tree branches.. We next indicate how to overcome this, problem. 

'-it somehow -* the voltages across r 4 and R8 controlled some dependent sources, 
they would become type 2 elements and would Hereford- 7b e branches of the 
















Figs. 2.13 and 2. lb: R1 in tree 

_ — - “ 








NASA P FROM LfSTgGj 


■ILTFR HOE 


NONE 

-NONE. 










loops. Conversely, the tree having the largest Branch Product (EL) yields the 
flowgraph with the greatest number of loops when compared with the flowgraphs 
generated by the R4 and R3 trees. However, the flowgraph generated by the R 1 
tree will give fewer loops than the flowgraph generated by the R2 tree since 
the Rl tree yields a smaller branch count. 

The circuit in Table 2.1 was analyzed using the mSAP program for the four 
possible trees available from the NASAP tree selection algorithm. The computer 
results are given in Figs. 2.9-2.16. : Fig. ^2 '.Ogives an input listing that 
includes r 4 in the tree. -Note that this is not the only listing order that 
will cause HI to be a branch of the tree --the only requirement on the input 
listing is that the card describing R 4 must precede those describing Rl, R2, 
and R3. Note that there are 332 flowgraph loops including 72 first order loops 
_ Fig.. >2.12 shows an input listing that makes R3 a tree branch and show that 
there are; 334 flowgraph loops including 74 first order loops. In Figv : ^2.l4' * 
is a listing with Rl in the tree . The flowgraph loops generated by the Rl 
tree total 366 including 80 first order loops ’ Finally^ Figv?* 2.16 gives a 
listing with R2 a tree branch. The flowgraph loops number 796 including l 86 
first order loops. Figs. 2.10, 2.11, 2.14, 2.15 give,', the transfer functions 
obtained for the different trees., 

Note that the R4 tree does indeed *yield the fewest flowgraph loops. How- 
ever there is almost no difference between the number of loops generated by the 
R3 tree and the R4 tree . The Rl tree yields about 10$ more loops than either 
the R4 tree or the R3 tree. The R2 tree generates more than twice as many 
loops as either of the other trees. 

A Needed Modification 

It should be noted that the primitive flowgraph as developed by the NAjSAP 
program does differ slightly from the actual primitive flowgraph developed 
from Kirchhoff’s voltage and current laws and Ohm* 


s Law. In the NASAP flow- 



graph there can be no connections to the current node of an independent or depen- 
dent voltage source as well as no connections to the voltage node of an independent 
or dependent current source. However such connections may exist in a true pri- 
mitive flowgraph. If there are no connections emanating from these nodes, then 
these nodes will not be a part of any loop and no information will be lost in 
determining the transfer function. ' This is the reason why one is not able to call 
for the current of a voltage source or the voltage of a current source as an 
output variable in the MSAP program. Furthermore, this difference between the 
MSAP and true primitive flowgraph affects the procedure for determining the 
optimum tree . 

Since there can be no connection between the voltage and current nodes of 
current and voltage sources (either dependent or independent), let us modify the 
branch count by omitting the branch count of those loops which are formed from 
independent and dependent current sources and by deleting from the branch count 
of each remaining loop those branches representing independent or dependent vol- 
tage sources . Similarly the Loop Count will be modified by omitting the loop 
count of those tree branches which represent independent or dependent voltage 
sources and by deleting those loops, which are formed by independent and depen- 
dent current sources from the loop count of the remaining tree branches. In 
other words, the modified branch count is the sum of the passive element tree 
branches that are part of those loops formed from the passive element links 
while the modified loop count is the sum of the loops, formed from passive ele- 
ment links, that pass through passive element tree branches. Note that the 
modified loop count will always equal the modified branch count. 

The need for the modified Branch Count, Branch Product, and Loop Product 
is demonstrated from the analysis of the Butterworth filter circuit in Fig. 

2.17.S [SA-1] i 
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Desired Transfer Function; W3/W1 

R3 R7 



Fig. 2.17: Equivalent Circuit of 

a Butterworth Filter 

The type 2 elements, V1-R2-V2-R6-V3 will be tree branches as will Cl and 
Ck due to the topology of the network. Either C2 or C3 but not both can be 
branches and the same situation holds for 05 and c6. There are four trees to 
choose from, namely 

1. with C2 and C5 

2. with C3 and C6 

3. with C2 and C 6 

4. vittx C3 and C5 

We will concern ourselves only with the first two choices. Table 2.7 
gives branch-loop data for the C3-C6 tree while the same data for the C2-C5 tree 
is shown In Table 2.8. 
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Table 2.7: Tree with C3 and C 6 

Link R1 R3 r4 C 2 R5 R7 r 8 C5 

32233223 

Branch VI V2 V3 Cl R2 C3 C4 R 6 C6 

143132132 

Branch Count = 20 
Branch Product = 1296 
Loop Product = 432 

Passive Link R1 R3 r4 C2 R5 R7 k 8 C5 

21122112 
Passive Branch Cl R2 C3 C4 R 6 C 6 

13 2 13 2 

Modified Branch Count = 12 
Modified Branch Product = 1 6 
Modified Loop Product = 3 6 

Table 2.8: Tree with C2 and C5 

Link R1 R3 R4 C3 R5 R7 R8 C6 

32233223 

Branch VI V2 V3 Cl R2 C2 c4 R6 C5 

132142142 

Branch Count - 20 
Branch Product - 1296 
Loop Product =384 

Passive Link R1 R3 r4 C3 R5 R7 R 8 c 6 

21222122 
Passive Branch Cl R2 C2 c4 R 6 C5 

14 2 14 2 

Modified Branch Count = l4 
Modified Branch Product = 64 
Modified Loop Product = 64 
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Figs. 2.18 and 2,19: Ti*ee with C2 and C5 
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tree. The problem remains of selecting these dependent sources such that they 
do not affect the circuit being analyzed. This is easily accomplished by making 
these sources dependent voltage sources with one node of the source connected to 
any node of the original circuit and with the other node left unconnected. Each 
unconnected node is actually an additional node in the modified circuit. Since 
a tree is an interconnection of all nodes and since the only connection to these 
hanging nodes is the source itself, then each of these sources must be a tree 
branch and thus a voltage source. However since one node of these voltage source 
is left floating, these 'sources will in no way affect the original circuit be- 
cause no current can flow through these sources. The zero-load feature of this 
type of "dummy" source is also quite apparent in the flowgraph. There is a con- 
nection to the voltage node of the voltage source from the voltage node of the 
controlling element. However no connection emanates from the voltage node of 
the voltage source 'since this "dummy" voltage source is not involved in any of 
the network loops formed by the links of the co-tree. A node can be part of a 
flowgraph loop only if there is a connection leading to and away from the node . 
Since these "dummy" voltage sources are not part of any flowgraph loop, then 
they cannot affect the determination of the transfer function. 

Thus the addition of the following two cards 

V4 1 11 1.0 mk 

V5 1 12 1.0 vr8 

to the input lists of either Eig. 2.18 or Fig. 2.20 will make r 4 and R8 tree 
branches without affecting the electrical properties of the original circuit. 
Nodes 11 and 12 are the floating nodes. Note that the dependency value (l.O' 
in this case) is completely arbitrary and can have any value except 0.0. In 
essence, the above "trick" simply is an artificial means to tell the NASAP 
program what elements we wish to have in the tree. It does have the drawback 
in that elements are needlessly wasted but this only becomes a factor when the 
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number of elements in the circuit is near the limit set by the length of the com- 
puter word. Note also that if an element which one desires to put into the tree 
by the above "trick” is also declared a type 6 element somewhere else in the in- 
put list (i.e., its current controls some source or is the desired output vari- 
able), then an error message will result. 

Fig. 2.22 gives the necessary branch data for the R4-R8 tree. Take note that 
the branch count, branch product, and loop product were determined without regard 
to the "dummy" sources V4 and V5. 

Fig. 2.22: Tree with r 4 and R8 

Branch count W l8 Modified Branch count, — ='12 

Branch product =576 Modified Branch product = 1 6 

Loop product = 2l6 Modified Loop product = 36 

Note that the modified data is identical to that obtained for the C3-C6 tree. 
However the unmodified branch count is two less than that of either the C2-C5 
tree or the C3-C6 tree. Figs. 2.24 and 2.25 give the computer results. The r4- 
R8 tree generates a flowgraph of 737 loops (including 17 first order loops). 

This is less than one -half the loops generated by the C2-C5 tree and more than 
200 loops less than the number of loops derived from the C3-C6 tree — a substantial 


reduction . 






CHAPTER III 


MODELING A CONTROL SYSTEM FOR NASAP 
IIIA GENERAL DISCUSSION OF CONTROL SYSTEMS 

The preliminary discussion of feedback control systems is kept "brief on the 
assumption that the reader already is generally acquainted with feedback control 
theory. Familiarity with introductory textbooks such as Dorf's Modern Control 
Systems [DO.-l^ and Perkins and Cruz Engineering of Dynamic Systems [PE-l] would 
be particularly helpful in that they use aerospace feedback control system pro- 
blems as illustrative examples. 

The major subdivisions of feedback control systems usually are: 

1. A plant, process or controlled system wherein the position or state 
is being regulated or monitored. 

2. The controller consisting of a sensor and control elements. 

3- A comparator or error- sensing, device to ^detect the- difference between 
the input reference and the output signal. 

Next we itemize the major steps involved in the design of a feedback control 
system. These steps are: 

a. Establish performance specifications for the system 

(e.g. type of control, tolerance on accuracy, speed of response, 
overshoot, etc . ) 

b. Interpret the specification data in terms of design parameters and 
components of the control system, (give due consideration to reliabilit 
space, cost, etc.) 

c. Formulation of the transfer functions of the components and analysis 
of the preliminary design. 

d. Improve the performance of the preliminary design by suitable com- 
pensation to meet the specifications. 

The scope of computer-aided design in this manual only covers aspects of steps c 
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and d. 



Control system design can be carried out either in the frequency or in the 
time domain. This is important in considering the possible role of NASAP in 
such applications. Control engineers have found it convenient, in the analysis 
of linear feedback control systems, to use the transfer function concept and the 
block diagram representation of the system. The transfer function concept is 
basic in the application of the frequency response method of analysis. In this 
approach the steady state response of the system to a sinusoidal input is used. 

The output/input signal relationship for each component of a control system is 
described by a transfer function. The operations of these components are repre- 
sented by noninteracting blocks which are interconnected to form the block dia- 
gram or the corresponding signal flow graph of the overall system. Thus one 
obtains a functional representation of the feedback control system equivalent 
to the set of simultaneous differential equations that relate the variables of 
the physical system. 

The basic procedures usually followed to analyze and design a feedback control 
system by the frequency response method are: 

i) Determine transfer functions for each of the components used in 
the system (from the differential equations via transforms or from 
physical measurements). 

ii) Formulate the signal flow graph from the system block diagram. 

iii) Reduce the complicated block diagram of the system to a simple single 
loop configuration having a transfer function for the forward and the 
feedback branch if the open loop transfer function and output trans- 
form of the control system are to be used. 

iv) Determine the system characteristics using the Bode plot or Nyquist 
diagram (an alternative graphical method uses the Nichols chart). 

v) To have the system meet the prescribed performance specifications, 
design the necessary compensators that will reshape the plots ob- 
tained in step iv. This may involve cut and try. 



In addition to the above, the designer often has to investigate sensitivity of 
parameters to variation of individual elements. The final step may include 
analog simulation or physical model tests. 

As an alternative to the frequency response approach, the analysis and syn- 
thesis of feedback control systems determination of the system stability and the 
evaluation of the output of the system in response to impulse, step; or ramp 
input functions. Here again there is often the need for compensation of the 
system so that it will meet specs. 

In the remaining chapters of this manual we shall indicate how HASAP can 
assist the design engineer to accomplish some of this work. It should be noted that 
some of these procedures are best carried out with a hybrid computer. Dr. C. H. 
Beck [BE-l] has developed a hybrid HASAP module for such applications as part of 
this cooperative development of the MASAP program. 

- _ t £* 

IIIB1 Equivalent Electrical Networks for Transfer Functions 

Before a control system can be analyzed using HASAP and before any necessary 
compensation can be determined, the dynamic characteristics of the plant must 
first be simulated by an electric network which has an equivalent dynamic charac- 
teristic. The transfer function of a lumped linear plant can be expressed as a 
ratio of two polynomials. The problem of modeling the plant transfer function 
can be simplified if the polynomials are put in factored form. The individual 
factors or group of factors can be modeled by using simple RLC circ ui ts . Then 
for a complicated transfer function these circuits are connected in cascade with 
suitable isolation between each circuit to prevent loading that would result in 
a change in the modeled transfer function. This necessary isolation is obtained 
by using ideal dependent voltage or current sources (which are available in 
HASAP) . Table 3.1 gives a list of some elementary circuits with isolation and 
their corresponding transfer functions. 
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IIXB2 Cascade Interconnection of Transfer Function Models 


As Indicated earlier the simple transfer function models can be cascaded. 
Consider a plant transfer function given as* 

, — ms+p ( 3 . 1 ) 

E^s) b(3s + l)(lts + O.ls + 1) 

This can be rewritten as: 


SgCs) 

i * 

( 2s + 1 \ 

/ 10 

) (3.2) 

SjW 

3 

^ 3s + 1 / 

' ‘‘S 2 + O.lB + 1 

Thus a ties cade interconnection 

of three circuits 

from Table 3.1 can be 

used to 


model the above transfer function. 

2 ,! 4H 



Note that Eg is a negative resistance and that the gain factor (for this example, 
10) is included In the dependency ralatlon for the third dependent voltage source 
(e^). The gain factor could Just as easily he included in the dependency relations 
of or e^ , 

Many of the subsequent examples used to Illustrate various aspects of computer- 
aided control system design will incorporate techniques for modeling the pertinent 
transfer function. 
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IIIC ADDITIONAL EQUIVALENT NETWORK MODELS 


IIIC1 Use of Negative R, L or C 

Since negative element values of E, L, and C are permitted In the NASAP 
input coding* rational transfer functions of control systems can he modeled 
simply hy application of the continued fraction expansion procedure. Accordingly 
the rational function will he represented in general by the input admittance or 
impedance of a ladder structure consistings of positive or negative R, L, and C 
elements. It is emphasized that this approach works because physical realizabi- 
lity as a passive network Is not a consideration. Only the equivalent dynamic 
characteristic matters. 

As an example, let us consider the biquadratic all -pass function 

r(s) = s ~ - ■ aa + b (3.3) 

s + as + b 

Performing the continued fraction expansion of F(s) yields 


s + as -t- b 


f 


s - as + b 


s 2 + 


+ t - aT fi 


-SaS 


as + b 


a -g 

as + b J -2as 


-2as -2b s 



~2b ) as + b 

as 2 

b Tab 


( 3 . 4 ) 


If F(s) Is assumed to be an admittance y , then the resulting ladder network 
is shown in Fig. 3*1* 
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In other words, the Input admittance of the circuit of Fig, 3-1 is a repre- 
sentation of the specified rational function F(s}* 

To utilize this Input admittance model with the NASAP program, two additional 
elements, a voltage source and a resistance of very small value, must be in- 
cluded in the circuit of Fig. 3.1 (see Fig. 3*2), 

l o 

2a Zb 



VIRl/VVI Fig , 3,2 

The voltage source VI represents the excitation for the input admittance 

Y- = ™ * The response is the current I , However the current flowing through 
V 1 1 

a voltage source cannot be specified with the TiASAP program. Since no element 
of the circuit of Fig. 3*1 is in series with the source VI, it is necessary to 
include the small resistor Rl* Thus the input admittance of this circuit can 
be specified as 

\ « na/wi (3.5) 

The circuit of Fig. 3*3 can be easily cascaded with other isolated circuits 
to model more complex rational functions. As an example, the transfer voltage 
ratio of the circuit of Fig, 3*3 models the function. 
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2 i 

3 - 3 + k 


(3.6) 


(s^ + s + U) (2s + 1) 



VV3/VV1 


S 2 - s + 4 
(s 2 + s + 4} (2s + l) 


Fig, 3,3 

Alternatively the function F(s) given above can be assumed to be an input 
impedance £ , From the continued fraction expans ion* we obtain the IttSAF 
applicable circuit of Fig, 3-4. 



Fig, 3*4 


The current source XI represents the excitation necessary for the input im- 
pedance, The response is the voltage across the very large resistance R1 which 
must be added to the circuit since no single element obtained from the continued 
fraction expansion is connected across the ideal current source terminals. 
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However it is possible to avoid the use of the small series resistor of 
Fig* 3-2 and the large shunt resistor of Fig* 3^ by taking the reciprocal of 
(i*e. inverting) the rational function that is to he modelled by NASAP before 
performing the continued fraction expansion* 

As an example, let us desire to model the biquadratic all -pass function given 
above by the input impedance of some ladder circuit* Ihus we desire a circuit whose 
input impedance is defined by 

(s) „ ^ - as * b ( 3 . 7 ) 

a + as + b 


Inverting this expression we obtain 


r a <0 


2 

s + as + b 
2 

s - as + b 


(3-8) 


Performing the continued fraction expansion this time yields the NASAP circuit 
is shown in Fig. 3*5, 



Fig. 3.5 


Note that now there is no need to Include a large shunt resistor across the II 
current source terminals (as in Fig* 3-*0 since the resistor R1 already shunts 
these terminals. 

Similarly, we can model the biquadratic all-pass function by the input 
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admittance of a ladder network defined by 


Y 


3 " 


2 

b - as + b 
5 

s + as + b 


(3-91 


Inverting this expression and then performing the continued fraction expansion 
(given above), we obtain the ladder circuit of Fig* 3*6, 



IRI/VVI = 


S z -os + b 
S z +os + b 


Fig- 3-6 


Note that addition of a small series resistor (as in Fig* 3*2) is not needed since 
the current flowing in the resistor HI is from the voltage source VI . 

In summary, by use of the continued fraction expansion procedure, we have 
obtained four NASAP-codable circuits to model the specified all-pass function 
F(s)* The circuits of Figs* 3*5 and 3.6 are more desirable as models than those 
of Figs* 3*2 and 3-^ since they require one less element* It should be further 
noted that each of these four circuits has another desirable feature* The IIASAP 
tree* selection algorithm will automatically select the shunt elements of each 
circuit as tree branches. This will lead to a star tree which generates the flow- 
graph with the fewest loops* 

This will always occur if the polynomials of the specified rational function 
are arranged in descending powers of s. The resulting continued fraction expansion 
will either make the series elements inductors or the shunt elements capacitors* 

T3 


IIIC2 Illustrative Examples 


Sometimes It Is necessary to arrange the polynomials In ascending powers to 
attain a ladder network with both positive and negative elements. For example, 
let us assume we wish to model the function 


V s ) 


~ 2 

s + 2s + s + 1 


(3.10) 


as the input impedance of a ladder network. We cannot perform a continued frac - 
— 1 1 o 

tion expansion cm (s) - s * 2s + s + 1, Consequently we will need a large 
shunt resistor in the NASA? model. Furthermore we must rearrange the denominator 
of F^(s) in ascending powers of s before performing the continued fraction expansion. 
The resulting NASAP model with shunting resistor R1 is shown in Fig. 3.7* Note 
that, since the shunt elements are 
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Z,- VRI/XJI- 


J 

S 3 + 2s*+S<H 


Fig. 3.7 


inductors, the HASAP tree -selection algorithm will pick a linear tree consisting 
of resistor R1 and three of the other four resistances (the particular R 1 s depending 
upon their location in the input list). This type of tree yields the largest 
number of loops in the corresponding primitive flowgraph. 

Sometimes in the continued fraction expansion procedure more than one term is 
eliminated by subtraction. This may necessitate a rearranging of the remaining 
polynomials to achieve a NASAP model. Such a case exists for the function 
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F 2 (s) 


+ 1 


( 3 * 11 ) 


The first three steps of the continued fraction expansion process is based on as- 
cending order of the denominator* Then three terms become zero, nils necessitates 
reversing the polynomial -s-s -s to -0 -3 -s. 

In Fig* 3-8 is given a MSAP-codable ladder network whose input impedance 
equals the desired function F^s), 


\ -\ l 



VRI/III = 

S 3 +S*+S + ! 


Fig, 3.8 


Another example Is the biquadratic function 


FjCs) 


2 

s + as + b 

^ 

s +■ as + c 


( 3 , 12 ) 


I*et us model this function as an input impedance* Inverting F^(s) and then per- 
forming the continued fraction expansion yields 


2 

s + as + b 


1 

n 

Is + as + c 

2 -JL, 

s + as + b c - b 

1 r 

(c - b) / b + as + s , 

/ C__- D 

b as 

Z - ] 

as + s / c - b 
c - b + 



( 3 - 13 ) 


7 



as + s 


T 
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c - b 


)• 


c - b 
a 

c ~ b 


The resulting network is shown in Fig* 3*9. 

b 

c-b 


-o 

e“b 



-q 

c-b 


( 3 - 13 ) 

cent. 


Z = VRI/IX I = 


S +o$ + b 


S + os + c 

Note the large shunt resistor is again avoided by use of the Initial inversion* Note 

also that in the underlined section of the continued fraction process it was necessary 

2 

to reverse the polynomial s + as + b since two terms were eliminated in the preceding 
subtraction. 


IHC 3 Equivalent Networks for Summing Element 

Since most control systems require some sort of feedback loop, an electrical net- 
work equivalent to the sumning (or subtracting) element that is compatible with NASAP 

must be used. Such networks are shown in Figs* 3* 10a and 3.10b* 

R 



Fig. 3.10a 

?6 


R 

■A/W 




V>- “(Vi + V!) 

Fig. 3 -10b 



III Fig* 3-lOa if <* = E si, the output voltage is equal to the difference of 
the two input voltage. Similarly If or = R = 1 in Fig. 1,1Gb, the output voltage is 
equal to the sura of the two input voltages. 
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II ID 


MODELS OF FEEDBACK CONTROL SYSTEMS 


IIID1 Examples of System Models 

Unis we now have all the elements necessary to model a feedback control system 
with an electric network that is compa table with NASAP« As an example, the unity 
feedback control system shown in Fig. 3.11* 


r(t) e (t) eft) 



has the following equivalent electric network 

R* 



Fig. 3-12 

where is equivalent to r(t) 

Vp is equivalent to e(t) 
is equivalent to c(t) 
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A second example shows a HASA? model for a non-unity feedback system (Fig- 3-13). 



Fig. 3-13 


Fig, 3-l4 is the equivalent electric network. 



Fig- 3-14 

where is equivalent to r(t) 

Vg is equivalent to eft) 

V _2 is equivalent to c(t) 

is equivalent to b(t) 

IIID2 Control System Model and Its Step Response 

Using a control system design problem adapted from D'Azzo and Houpis, pp, hOS- 
4ll LM-l] we illustrate modeling of the system and use of NASA? to tabulate and plot 
the step response. 

Consider the unity feedback control system with cascade lead compensation shown 
In Fig- 3*15- To determine the step response when K = 10 we first obtain the MASAP 
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f i ' 

% 

Fig. 3-15 


circuit model show in Fig- 3-l6- In this model 
V5 represents the input R(s) 

V4 represents the output C(s) 

V6 represents 2.06 times the error signal; 2.0 6 E(fi) 

Ifote that the system gain* K = 10* is included In the dependency relation of VI* 
Also In Fig* 3'l€ we have NASA? output* namely the transfer function zeros and 
poles based on a simple flowgraph having only right first order loops* The de- 
sired step response is given in Fig* 3-17- Figure 3-I7& Is the response func- 
tion and 3.17b gives 51 discrete time values from zero to six seconds in equal In- 
crements. Finally in Fig* 3-18 ve have the plot of this time response* 

These NASAP outputs can be used to determine the key step response charac- 
teristics, Note in particular that the steady state value is unity since 
Fijff. 3*l?a we see that the residue of the pole at the origin is unity. Further- 
more in Fig* 3*17b the peak overshoot occurs at t « 1,^ seconds and is 9,15$. 

This response settles (for 2$ tolerance) in 2,04 seconds. 
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NASAF MODEL KB CONTROL SYSTEM 





Fig* 3.16 

Resistance Values in Ohms 
V5 represents the input R(a) 

V4 represents the output C(s) 

V6 represents 2*06 tijnes the error signal; 2,0 6 E(s) 
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CHAPTER XV 


CONTROL SYSTEM ANALYSIS IN THE FREQUENCY DOMAIN 
IVA ANALYSIS OBJECTIVES 

As mentioned in Chapter III control system design is usually aimed at a set 
of specifications. Therefore we start with a summary of typical performance spe- 
cifications as general "background. For convenience we adapt the tabulations 
given by Grabbe, Ramo, and ^Woodruff :[GA £] and list some of them for transient 
response in Table 4], la and for frequency response in Table 4.1b. These are suf- 
ficient to indicate the sort of specifications that can be expected for the class 
of linear time -invariant single input single output control systems being con- 
sidered in this manual. 

There are several approximations or rules of thumb which were developed by 
control engineers for use ‘.when time or facilities are not available for a more 
exact analysis .of this class of systems , They are also useful as rough checks on 
the results of a computer analysis. The more common of these rules of thumb are 
presented, virtually unchanged from [GAS] in Table 4.2. ^They'mustj, however, 
be used with caution since being approximations , they do not apply with equal val: 
dity to all control systems. Note that the approximations for transient response 
are applicable only for step inputs. 

Following Table 4.2 we have Fig. 4.1 which shows typical step and frequency 
response curves to help pin down the definition of some of the terms used in 
these tables. 


85 



Table 4.1a 


Common Performance Specifications 


Type 

1. Transient oyershoot 


Definition 

Usually taken as ratio of 
peak of transient to final 
value for a step -lop at. 


2. Settling time Defined as time to reach 

and remain within a spec- 
ified percentage of final 
. value (often as 5% or 
2%) after a step input, 

3. Steady-state error Final error existing between 

desned and actual out- 
put. 


4. Rise time 


- *'5. Dead time 


G. Absplute damping, 
decrement factor 


Defined as (a) time to \ the 
final value, or (6) slope at 
\ the final value, or (c) 
time between 10% and 
90% of final value after a 
step Jnput, 

* Defined as (a) time for out- 
put action to be initiated, 
or (&) f6r output to reach 
a given level (10% or 
50%), or (c) -time to the 
intersection of the sldpe 

of the transient at \ the 
final value and the initial 
value after a step input. 
; Defined as the real pait of 
' the roots of a quadratic 

system and as such deter- 
mines the rate of decay of 
1 bans ient. 


7. Damping ratio » Damping ratio is defined as 

f in the quadratic s 2 -{- 
2&? 0 fs H- too 2 and indicates 
the decay per cycle of the 
natural fiequency. 


Transient response 

General Remarks 

Convenient when transient solution is avail- 
able. Can be estimated fiom loot locus or 
frequency response. Useful foi nonlinear 
’systems System must be j ~ excited by 
^ s tep inputs and be undei.dam.ped J " J 

See 1 above Used for systems which i equirc 
rapid synchronization, eg., file conti ol sys- 
tem. 


See 1 above. Easily calculated fiom static 
characteiistics or final value theoiem. Use- 
ful when input is simple aperiodic function. 
Can include frequency components which 
arise in nonlinear systems. 

Easily estimated fiom frequency response or 
root locus and is indicative of band pass of 
system. Used for overdamped systems. Has 
found application in process controls where 
characteristics (1) or (2) may not be easily 
recognized. 

See 1 above. Easily estimated from fiequency 
response and is indicative of phase shift near 
gain ciossover in systems. Useful when de- 
lay times exist in system. Used for over- 
damped systems. Both rise and delay time 
derive from filter theory. 


Convenient method of interpreting more com- 
plex systems in terms of quadratic systems. 
Valuable in combination with relative damp- 
ing in work with root locus analysis. Has 
had extensive use in systems demanding pre- 
scribed transient peiformance, paiticularly 
when the time decay is important, e.g., in 
autopilots. 

Useful because it is a parameter in nondimen- 
sional plot of quadratic response. Used in 
combination with 6 above in root locus anal- 
ysis. Used when number and size of pver- 
shoot are important In combination with 6 
above defines decay of oscillatory component 
of transient. 


86 



Table h.lb 


| Common Performance Specifications ■ Frequency response- 

8. Phase margin DefineHas lSO 0 + phase shiT b Used as" a rule of thumb in frequency response 

at unity gain of the open analysis to indicate stability and peifoim- 

* loop f i equency response. ance. Easy to use and to obtain dii ectly from 

• frequency lesponse diagi am. 

9.. .Gain margin Gain maigin is ratio of max- Same as 81 Indicates relative sensitivity of 

imum stable gain to ac- system to gain vaiiations. Can be calcu- 

tual gam, i e., gam at latcd by Pouth’s criterion. Not as good a 

, * phase crossover. criterion for performance as S. Little used, 

peak Patio of maximum of closed Used with Nyquist and frequency response 

loop frequency response analysis. Rules of thumb i elate and 

to a low frequency value. transient ovei shoot. Easy to calculate from 

frequency lesponse diagram. 

If. Band width Defined vauously (a) usu- Used i\ith fi equency lesponse analysis and is 

. 1 ally as frequency where related to speed of response of system. Used 

closed loop response falls also when definite frequency bandpass is 

to or 3 db of its low f re- needed for fidelity. M m , bandpass, and the 

quency value,or (6) some- phase shift at these valves give a good indi- 

tim'es as the frequency at cation of the closed loop response and are 

the significant peak M m , often used when a number of closed loops are 

or (c) the crossover of the operated in tandem as system, 

open loop response. 

12* Static error coefficient "i Defined as the final eri or re- Used to set low-fi equency gain of open loop 
suiting from a continu- frequency lesponse. Useful where steady in- 
ous input of position, or puts are encountered, 
velocity, or acceleration, 
etc. The magnitude of 
the input and the maxi- 
mum toleiable enoi must 
be specified. 

13. Dynamic error coeffi- Defined as the steady-state Relates system gain and time constants to er- 

cients (or steady-state eu or resulting fiom the rors aiising from liighei deiivatives of input, 

enor coefficients) derivatives of the input Used to estimate erroi lesulting from varying 

function. The time func- input to given system and conversely to de- 
tion and/or its deriva- tci mine closed loops pole-zero location to give 

tives must be specified as desired error. Accuiate wheie input varies 

well as the maximum tol- at slow rate compai ed to bandpass. Becomes 

erablc error. pooler as input vaiies more lapidjy because 

of tiansient effects. 
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Table k.2 


Paiameter 
Time to 
peak 


Peak over- 
shoot 


HuLE-OF-TlIUiU B APPROXIMATIONS 


Approximation 
l P ~ tt/wc 

where t p = time from step in- 
put to peak value 
of response tran- 
sient, seconds 
w e = open loop crossover 
frequency, ladians/ 
second 


C/R\ p ~ 0.85 M m 
wheie C/R\ P ~ peak value of 
transient re- 
sponse to a 
step input 

M m — maximum val- 
ue of closed 
loop fiequency 
response 


Remarks 

InC4£,e5nvith a dominant complex 
pair of closed loop poles, the 
open loop ci ossovci fiequency, 
w e , times the time to peak, t P) 
is about 3 or ir. In other woi ds, 
the time to peak is about half 
the period corresponding to the 
^open I ooj^ crossover frequency. 

The peak value of the transient 
response, C/R\ p , to a unit step 
input is geneially less than the 
maximum steady-state value, 
M m , of the closed loop fie- 
quency response The maxi- 
mum value of C/R\ p generally 
approaches 2 0 while the maxi- 
mum value of. M m approach- 
es infinity. For many applica- 
tions "good” servos are those 
with the values of M m be- 
tween 1.3 and 1.5. 


5 Damping 
ratio 


Settling 

time 


r = i/C2 m c ) • 

where f = damping ratio 
M* = value of closed loop 
fiequency response 
at the corner fre- 
* quency 


6%) « 3 VF- 
« 5Vl - r 2 /fw d 

where t, — time for response to 
step input to settle 
to within some per 
cent of final value, 
seconds 

Teq = time for response to 
reach 63% of final 
value 

o>d - damped natural 
frequency, radians/ 
second 

f = damping ratio 


The damping ratio may be ap- 
proximated from the value of 
the closed loop frequency le- 
sponse of the system at the 
comer fiequency, (the fre- 
quency at which the lines 
asymptotic to the log magni- 
tude cuive intersect). This is 
exact for a’second order system. 

The settling time, t a , is generally 
defined as the time for the sys- 
tem to settle to within 5 or 

sometimes 2% of the final val- 
ue. In either case it is quite 
difficult to predict t, for an 
underdamped system because 

. it is subject to fluctuations of 
about one-half the period of 
oscillation for onty small 
change s in system parameters. 
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Time response 


Table 4.2 


Rule-Of-Thumb Approximations ( Continued ) 


Parameter Approximation 

Rise tim e i^t ~ tr&m ~ 1*3 

where f r = rise time (10 to 
' - 90%) 

a t = (dtefined above) 
ca m = (defined above) 


Remarks 

The system’s rise time, t r , which 
is here considered to be the 
time for the response to a step 
input to go from 10 to 90% of 
its final value may be approxi- 
mated as indicated for systems 
with a Afm value of about 
1 3 to 1.5. 


Phase mar- 7« = 40° 

gin at where 7c = open loop phase 

crossover margin at the 

frequency crossover fre- 

quency 


A phase margin of 40° at the 
unity gain (crossover) frequency 
generally corresponds to a 
M m ratio of approximately 
1.5. Since this value of M m 
is the maximum ordinarily con- 
sidered feasible, the phase mar- 
gin should be 40° or gi eater. 


Oscillation 

frequency 


0}( ~ CJ m ~ 0.7O6>e 
where at — oscillation fre- 
quency of tiansient 
response, radians/ 
second 

(o m — frequency at which 

M m occurs, radi- 
ans/second 
a c = open loop gain 
crossover fre- 
quency, radians/ 
second 


The frequency of oscillation of the 
transient response, gu, is gener- 
ally about equal to the fre- 
quency, dm, at which the fre- 
quency response peak, M m , 
occurs. Both a m and are 
usually less than a e , the open 
loop crossover fiequency. For 
the “good” servos with M m 
= 1.3 to 1.5 an approximate 
relationship is as indicated. In 
this approximation at is used to 
mean essentially the same thing 
as a d , the damped natural fre- 
quency, previously defined for 
a system with a dominant com-: 
plex pair of poles. J 




Figure, 4 J-- (a) ^Typical serve system response to unit' step input. 

(k) Typical system frequency response (close d.loo p). 
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B7B BODE AM) ROOT LOCUS PLOTS 

Bode showed that the phase angle of a (minimum phase) transfer function G 
could be related to the rate at which the magnitude of G decreases with in- 
creasing frequency. This is the basis for the "frequency response" j, method of 
analysis. In this method, the magnitude of G in decibels and the phase angle, 
are plotted on semilog paper as functions of the frequency a) (plotted on the log 
scale) with jo) substituted for s in G. The value of |g| in decibels (dB) can be 
found from the value of |g| by the following equation. 

l G ldB " 20 log 10 l°l (*►.!) 

For preliminary estimation of compensating networks to insure stability or im- 
prove performance;, it is often possible to omit the phase angle plot and to 
make use of an approximate attenuation plot (or plot of |g|). Exact plots 
are needed for a final check after selection of the proposed stabilizing transfer 
functions, however. The method of drawing approximate attenuation plots is des- 
cribed in most textbooks on control theory. Before we discuss and illustrate 
how NASAP obtains exact plots of the magnitude and phase of G, we provide a col- 
lection of typical transfer functions with their corresponding Bode plots in 
Table 4.3 » These plots indicate the gain and phase margins. 

Another important tool for analysis and synthesis of linear control systems, 
usually attributed to Evans, is know as the "Boot Locus" method. As with the 
frequency response methods, its importance derives from the fact that it helps 
provide insight into the significant aspects of any particular system. It is 
not restricted to direct feedback systems nor to systems with open loop poles 
and zeros in the left half plane. 

it is recalled that for any closed- loop system with "degenerative feedback" 
see Fig. 4.2, the closed-loop transfer function is 



c 

R 


G 

1 + GH 


(k.2) 



Fig, 4.5: Block diagram of a basic 

feedback control system. 

The closed-loop response is determined by the roots of the denominator, i*e., 

the characteric equation. The roots must all be in the left half plane in order 

that the system be stable. Furthermore the time-domain design parameters, such 
as peak time maximum overshoot, damping factor, and settling time, are initi- 
mately related to the s -plane location of the roots of the characteristic 

equation of the control system, which are the poles of the closed- loop system 

function , Accordingly, a knowledge of how the roots of 1 + gH vary when the 
gain constant of GH varies should be of considerable assistance in understanding 
the system. A plot of the locus of the characteristic roots of the control sys- 
tem with the system loop gain as a parameter is commonly known as the root locus. 
The design of feedback control systems by use of the root- locus method 
involves the reshaping of the root- locus plots by shifting or introducing open- 
loops poles and zeros* As a preliminary to the discusssion of the design 
aspects, the effects of shifting open-loop poles and zeros are first indicated 
with the aid of examples in Table 4,3. 

Gain Adjustment ■ When the preliminary analysis of a control system indicates 
that the system is unstable or that the over -all performance is inadequate, 
steps must be taken to improve the system performance. The most direct and 
simplest way of changing the performance is by the adjustment of the system gain. 
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Table 4.3 


Plots For Typical Transfer Functions 
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(Continued) 
















However, for most control systems the design specifications cannot be met by 
gain adjustment alone. The usual alternative' is the introduction of compen- 
sating devices into the control system. 

The adequate gain setting for a control system can be determined from the 
gain-phase plot of the system or from the Bode diagram of the system. 

~ change in system 'gain usually affects practically all of the system de- 

sign parameters.. For instance, an increase in system gain may cause a reduction 
of the system error, may increase the speed of response of the control system, 
and may make the system more oscillatory. The effects of gain variations upon, 
the behavior of a control system are conveniently observed on the root locus 
plot of the control system -also. 

NASAP Output 

For a given transfer function G with all the coefficients known, the fre- 
quency response is calculated by setting s « jto and simplifying the expression 
to a linear combination of real and imaginary terms: - 

G = A (to) + jB(co) . (4.3) 

The magnitude of G and the angle 8 (to) are computed according to the equations. 


|G(jto) j = 

A 2 (to) + B 2 (to) , 

(4.4) 

9 (to) = tan 1 

B(to) 

AW * 

(4.5) 


Now if to is made to vary, then for each value of to the |G(jto)| and 0(to) can be ^ 
obtained ,‘bver the frequency range of interest and thus can be made ‘.available for 1 
plotting. With the complex arithmetic capability of FORTRAN IV, these computations 
are easily done in NASAP. The Bode plot consists of the |G(jto)| in decibel units 
and 0(o>) in degrees versus the log^to, taken over the frequency range specified 
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■by the user. 

To obtain a root locus plot for a control system, Fig. 4.2, it is required 
to find the values of s for which GH = -1 (or 1 + GH =0). For this it is ne- 
cessary that the angle of the complex number, GH, be 180 degrees and the magnitude 
of GH be unity. Thus, the complex number, s, must be selected so that the angle 
of the complex number, GH, is l80 degrees. When such a complex number for s is 
determined, a value of gain K can then be found which will make the magnitude of 
GH unity, although this value of JC might not' necessarily be 'the ?same, as the value 
specified in the transfer function. However, after a locus of values of s for 
which GH = 180 degrees has been found, somewhere along this locus one can find 
a number that yields |gh| = 1 f or the specified value of K. MSAP furnishes the 
necessary data in tabulated form to obtain the locus of points for which GH = l80 
degrees. To obtain the root locus plot directly by the computer requires an extra 
program that mil not be described here. Such programs are available in the lit- 
erature as illustrated by Program D91RTL by Vernon [Appendix I in VE 1] and 
another by Krall and Fornaro [see K& 1 or KR2]. 

An alternative approach is to use the root sensitivity data available from 
ITASAP to approximate the root locus plot. Such sensitivity data for an aerospace 
control problem is given in Chapter VI. 

As an illustration of the open loop Bode plot output of NASAP we use a con- 
trol system plant whose transfer function has a seventh degree polynomial denomina- 
tor. The NASAP print out is shown in Fig. 4.3. This Eisenberg control problem 
[El 1] is discussed further in Chapters V and VII. 
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NAS AP ”P~R OBI’S M E’fifE N BE R G “C'Q N T tfO l SYSTEM 
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IVC USE OF COMPENSATION 


When the desired behavior of a control system cannot he obtained by the 
gain adjustment alone, compensation techniques must be used. Compensation 
means to improve the system performance by reshaping the open-loop transfer 
function characteristic of the system. The compensation of a control system 
can often be accomplished either by an element in series with other components 
as shown in Fig. 4.4a or by an element in parallel .withone or more components' 


Compensating 
| networks j 

i 1 

G c {s) 


Controlled 

system 


G^s) 


C(s) 


Fig.~~4.4a Cascade compensation 


to form a subsidiary loop, as shown in Fig.- 4.4b. The former arrangement is 
referred to as cascade or series compensation and the latter is called feed- 
back or minor -loop compensation. A compensator or compensating device can 



Fig. 4.4b Feedback compensation 


stabilize a system which is unstable for all values of gain; it can improve 
both the transient and the steady- state performance of a system; and it can 
reduce the system error. Compensators are classified according to their 
operating characteristics into phase -lead (differentiating) type, phase-lag 
(integrating) type, and lag-lead (integro-differentiating) type. An example 
of the latter was given at the end of Chapter 3. The phase-lead type of com- 
pensator is generally used to modify the high-frequency portion of the open- 
loop transfer function plot and to improve the transient behavior of the system. 
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whereas the phase-lag type of compensator is often used to alter the low-fre- 
quency portion of the open-loop plot and to improve the steady- state performance 
of the system. The lag-lead compensation gives results intermediate "between 
these two extremes. By proper adjustment of time constants, considerable 
flexibility of control system characteristics may be obtained. 

Th choice of a method of compensation generally depends upon the specific 
system involved, the available components, cost considerations and the designers’s 
experience and judgment. 


The most commonly used, configuration for compensation is shorn in Fig. 4.4a 
The procedure to compute G Q (s) is to first compute the open loop transfer func- 
tion and then compute G c (s). There are two disadvantages in using this confi- 
guration. First if the overall open loop transfer function G^(s) is not pro- 
perly chosen, the compensator computer in this manner may not be realizable as 
a RC network. Second, it generally requires pole-zero concellation. In sys- 
tem theory terminology, it means that some poles of the overall system are un- 
controllable and/or unobservable. This poses a design problem in that these 
poles are dictated by the given plant G^s) and cannot be controlled by the 
designer. 


The feedback compensation 'y jshown in Fig. 


4‘.-4b, is sometimes superior to- 

_ " n 


cascade compensation in that variation of the parameters of the system compo- 


nents bridged by the feedback elements of the minor loop have less effect upon 
system performance if the minor-loop gain is made sufficiently large and if the 
parameters of the feedback compensator do not vary. Although the • 
freedom in choosing (L (s) and H (s) is greatly increased, the difficulties en- . 
countered in Fig. 4.4a may still occur. Therefore the series and shunt com- 
pensation are not always applicable in practice. 
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Consider the configuration of compensators introduced hy Chen [CH l]. 


This is shown in 


W 1 (s) 


Fig. 4.5* where k is a constant gain, C^(s) - D 


and 



Fig. 4.5 Chen’s Control System Configuration 




C 2 (s) = D ^ are proper rational functions with the dgree of polynomials 
I>l(s), D 2 (s) and N 2 (s) equal to n-1 and that of H (s) - n-2. This control 
system configuration has the overall transfer function 


kG(s) 

= 1 + c-^s) + kc 2 (s; g(s; 


(4.6) 


aro 2 d x 

= D D 1 D 2 + W 1 D 2 D + Mf 2 H D x 

If the denominators of C^(s) and C 2 (s) are chosen to he the same, that is 
Di(s) =Dg(s), then the last equation reduces to 


m>i 

G f ^ “ DD~ ""+ l^D + kN 2 N 

It can he shown that hy using this system configuration the compensators can 
he always chosen to he realizable hy RC networks and the cancelled poles can 
also he controlled hy the designer [see CH l]. It is worth noting that the 
complexity of compensators in the system of Fig. 4.5 is comparable to that 
required for the compensator of the corresponding control system based on the 
configuration Of Fig. 4.4a. 
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Examples 


An example of cascade compensation -was included at the end of Chapter III 
(refer to Figs. , 3^16 and 3 * 17 )- Here we provide the Bode plots in Figs. 4.6 and 
4.7 for the uncompensated as well as the compensated case respectively 'jbo faci- 
litate comparisons. 

Next we illustrate basic feedback compensation with the aid of two exam- 
ples having the configuration shown in Fig. 4 . 8 . The feedback function H c is 
used to modify the characteristics of the plant G, while the cascade function G c 
'-Is provided to aid in adjusting the performance of the major loop. Often the 
cascade compensation G c is a simple gain factor used to adjust the degree of 
stability of the system. The burden of modifying the transfer function G is 
placed on the feedback compensation H c . In addition, the feedback function is 
usually provided with an adjustable gain factor to permit setting the degree of 
stability of the minor loop. 



Fig. 4.8 General arrangement for feedback compensation. 

The procedure for adjusting the feedback compensation can be based on the 
magnitude asymptotes of the minor loop. In the frequency ranges where the open- 
minor-loop frequency-response magnitude |G(jo>) H c (jo>)| is very large, the closed- 
minor-loop response |c(jco)/a(jco) | behaves like the reciprocal of ^the feedback 
compensation. When the open-minor -loop response magnitude is very small, the 
closed-minor-loop response behaves like the plant. 

The feedback compensation is jused primarily to improve the dynamic behavior 
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of the plant in the mid- frequency range without altering the high gain of the 
plant at low frequencies . A high-frequency boundary will always exist since 
in any practical case the magnitude of G(jco) H c (jco) will become less than 
unity as frequency increase. 

The procedure for adjusting the feedback function H c and the cascade gain 
factor [G c (s) = K c ] can thus be roughed out by means of asymptotic plots of the 
pertinent responses and then carried out in detail by means of the gain-phase 
or Bode plots. The asymptotic plots enable one to examine the form of the 
closed-minor-loop response as the feedback compensation is adjusted. 

_ 5 Since there are usually several' parameters to adjust in the feedback 
compensation procedure, the process of design is one of trial and error. 
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For the above unity feedback control system with lead- _ compensation, get Bode- plot.- 










This example, shown in Fig. 4.9 along with its NASAP circuit model, is 
taken from Newton Gould and Kaiser [ME 1 pp. 326-333]. One finds that the 
magnitude plot of the closed -loop transfer function shows no sign of resonance 

I 

effects so that a reasonable closed-loop performance may be expected. The 
corresponding phase curve is given in Fig. 4.10. 
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Fig. 4.10 Closed-minor-loop phase-angle response. 

Using the 45° phase -margin criterion to adjust the cascade compensation, 
the magnitude crossover frequency of the major loop is found to be 

= 7*7 rad sec" 1 (4.8) 

and the corresponding gain factor of the cascade compensation is 

K = 38 ( 4 . 9 ) 

If this control system exhibits undesirable performance because of the lag-com- 
pensation effect in the open-major-loop response at low frequencies, some 
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Improvement can be expected by decreasing the feedback compensation time constant 
T c and by some increase in the minor-loop gain factor K q . The degree of improve- 
ment achievable must be ascertained by further UASAP runs . Further discussion 
-of this example is found in Chapter V including step response and error response 
for the control system with its gain K = 38. 
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0.16499976 01 0*44668036 02 0,40313256-01 -0.818O165E 00 0, 1004652E 01 0.20156626-02 

l CU46Q3Q3-7E-Q-1 0-564-1-8386— 02 0, 2496322E=01 ^a.5235.133£-00 4.4QQ28J-&-E— 0.] 0.424348.6 E=02- 

0.17499976 01 0.56233766 02 0. ) 5517706-01 -0.4100884E 00 0.10017B8H 01 0 . 7758053E-O3 

i 0 ^. 1 . 79099 . 70 - 0-1 0nr6 3095356-02 0 .-967-8025 E-02 =4. 230144 J-E-OQ tU4aQ4U-5£-UJ 0 .-43390 12E=03- 

] 0.18499976 01 0.70794026 02 0 .6044853E-02 -0.20527896 00 0. 1000696b 01 0 . 3022428F-03 

Q.18.990Q.7E .01 Q.7943224E 02 Q. 37764576^02 *-0.145^4166 QO 0.100Q4 35E 01 0,4 38823.93-03- 

0.19499976 01 0.89124476 02 0.23687606-02 -O.1027696E 00 0. 10002736 01 0. 1184380E-03 

l_U 0-r*l-999997 E-0-1 q^99-92934 6—0-2 O-.440 262-1 H^O 2 0 .-7.2423-673^01 0^1,000441 E-U4 0,44-13 1,07 6»0A- 








The configuration of the second multiloop example is also similar to that 
shown in Fig* 4.8. In this case however the G c block is simply* a direct con- 
nection as shown in Fig. 4.11. Figures 4.11a through 4. Ilf include the NASAP 
printout of the Bode plots and the step response. 
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CHAPTER, V 


CONTROL SYSTEMS ANALYSIS IN THE TIME DOMAIN 
VA INPUT SIGNALS FOR TIME RESPONSE' ' ~ 

The impulse function is the basic tool for analysis and synthesis of linear 
systems. However specifically for the study of linear control systems the 
unit step function response is the most widely used with the unit impulse and 
the unit ramp functions also commonly used as test inputs. 

We- note in passing that these three functions are related to each other by 
one or more integrations or differentiations. For example the unit ramp as a 
function of time is the integral of the unit step function. On the other hand 
the unit impulse may be regarded as the derivative of the unit step function 
(this concept is adequate for the purposes of this manual even though this 
derivative does not exist in the sense of calculus). 

One of the first steps in adapting NASAP to control system design was to 
investigate the feasibility of incorporating additional input functions. This 
is discussed next. 
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VBv',' A gPlT-IQML CTFU^S FOR CONTROL APPLICATIONS . 

Although the NASAP 69 /X package provides the time response of the out- 
put only for an impulse excitation, the time response subroutine INV of 
NASAP can easily be extended to the other basic input signals often used 
in control theory. The necessary additions to NASAP incorporate the 
step and ramp excitations are given in Appendix A. 

The algorithm used by NASAP to determine the residues of the poles of 
the transfer function assumes that the poles are simple. This poses no 
problem since the root-finding subroutine of NASAP will not find double or 
higher roots but will locate a number of simple roots iruthe neighborhood 
of the actual higher order root location. Thus for an actual double root on 
the negative a-axis, say at 

S l,2 = - a > 

the root finding algorithm of NASAP will indicate complex root pair at 

S 1 2 = a ± 

where the imaginary parts of 2 are extremely small in comparison to the 
real part. Thus while the analytic time response for a system with higher 
order poles will not be correct, strictly speaking, the table and plot of 
discrete time values will be sufficiently accurate for practical purposes. 
Given a rational transfer function 

'’"h(s)=§M (5..1) 

and the Laplace transform' of the excitation R(s) , the Laplace transform of 
the output can be expressed 

C(s) 53 H(s)R(s^2 

Then, by finding the residues of the poles of H(s)R(s), the time response 
of the output c(t) is readily obtained. 

If the excitation is an impulse function, R(s) =1, then the Laplace 
transform of the output C(s) is numerically equal to the transfer function 
H(s). This is method used by NASAP to find the impulse response. Once the 
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transfer function and its poles have been determined, all that is necessary 
to find the corresponding time response of the output is to calculate the 
residues of the poles of the transfer function. 

However, if the excitation is a unit step, R(s) = , then the poles 

of the Laplace transform of the output are the poles of transfer function 
H(s) in addition to a pole at the origin. If H(s) has no higher order poles 
and no pole at the origin, then the analytic expression and the tabular 
values of the output will be correct. Note that the residue of the pole 
at the origin is the steady-state value of the output. If H(s) has higher 
order poles but no pole at the origin, then even though the analytic expressic 
for the output as determined by NASAP will be somewhat in error, the tabular 
values will be correct. If H(s) does have a pole at the origin, then the 
results obtained from NASAP for a step response will be in error and will - 
probably cause premature termination of the execution due to division by zero. 
This can be seen as follows : 


K(s+a)(s+a 9 ) ... (s+a) 
J. Z — m 

s (s+b- ) (s+b 0 ) ... (s+b ) 
1 Z n 


Then 

C(s) 


^ K (s+a- ) (s+a ) . . . (s+a ) 
1 Z m_ 

s s(s+b-) (s+b 0 ) . . . (s+b ) 
1 Z n 


when R(s) = — . 


( 5 . 4 ) 
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residue approaches infinity. This situation can he avoided by the addition 
of some elements to the original circuit to put a zero at the origin and there- 
by determining the ramp response of this new network. In Fig. 5*1 is given a 
two port network W with a source V (either an independent current or voltage 
source) and an element x whose voltage is the required output quantity. The 
transfer function of this network has a pole at the origin and it is required 
to find the response -of the voltage across the element x to -a unit step 
excitation . 



Fig. 5.1 


Wow consider the modification of the network in Fig. 5.1 shorn in Fig. 5.2. 



Fig. 5.2 


where is a dependent voltage source whose voltage equals that across the 
element x and is a one farad capacitor. Wow 



V s ) 

sC aV s) 

sV A (s) 



~WJ 

" V(s) 

_ V(s) ' 


But 

V s > 

= V (s) 

X V 7 



so that 

I A (s) 

sV x (s) 

sP(s) 

P(b) 


TtJJ 

' V(s) ' “ 

sQ.(s) 

ottt 


Thus the pole at the origin is effectively eliminated. The response of 
i (t) to a ramp input is equivalent to the response of v x (t) to a step input. 
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A similar method can be used if the output quantity is a current through 

an element x. In this case is a dependent voltage source whose voltage 

equals the current through the element x. 

If the excitation is a ramp, r(t) = t, then R(s) = ~~ • Thus the 

s 

poles of the Laplace transform of the output are the poles of H(s) in addition 
to the double pole at the origin. As noted earlier the algorithms of 
NASAP can only approximate double poles. To illustrate this assume that the 
double pole at the origin is approximated by a complex pole pair located 
on the j axis a distance) aifrom the origin. Thus 


R(s) (s+j a)(s-g a) 


(5.8) 


partial fraction expansion yields 


R(s) 


1 1 
2j a s-j a 


i_~ 

s+j a 


Taking the inverse Laplace transform yields 


(5.9) 


. v 1 , i at -lat. 1 

r (t) = ttt— (e J -e J ) = — sm at 
a a 


air e ' 


j at -j at 


+ j 


2 a 


(5.10) 


However if at is .'.small, sin at can be approximated by- at. This 

... 1 


i(t) 00 t 


if 1 at « 1 


Thus the smaller) ra. : is, the larger the range of time values for which the 
approximation closely resembles a ramp input. 

Note that the coefficients of the exponential terms in (5.10) are 
conjugate imaginary. However, in general, these coefficients will be com- 
plex conjugate. The significance of this is now demonstrated. Let us 

i a t “j ct‘‘t 

assume that the coefficient of e J is A-jB while that of e J ' is A+jB. 
Thus we have 
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CA-W at + (A+jB) e - ^ 0,t = ACe 3 0 V e _;i at )-jB(e^ at - e" j a e j 

- 2A cos a^t + 2B siti’out (5.11) 

Now let 

2k = M sin 6 '(5.12) 

2B - M cos 0* (5.13) 

Substituting (5.12) and (5.13) into (5.11) and then simplifying yields 

(A-jB)e^ a Z + (A+jB) e 5 at = ^sin < at + 8) (5.14) 

where 

M = 2\! A 2 + B 2 ‘ (5.15) 

and 

0 = tan -1 ~r (5.16) 

£> 

Now if A « B, then (5.15) and (5.16) can be approximated by 


M ^ 2B (5.17) 

e- » f- (5.18) 


Since 0 and< a- are both small, there is a range of values of t where 
the approximation 

sin( at + 0) at + 0 

is valid . Thus 

(A-jB)e j at + (A+jB)e~ j at ^ 2 tx Bt + 2A (5.19) 

for A << B which represents a ramp with slope 2 aB plus a step function 
of magnitude 2A. 

If the constant A is zero, then the above equations approximates only 
a ramp of slope 2 aB. This agrees with equation (5.10). On the other 
hand, if A » B, then (5.15) and (5.16) can be approximated by 


M 2A 
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Thus sin)( at+ <f> ) ^ sini( at + y ) cos at. Since} at is quite small for 

a wide range of time values, then 

cos at#: 


This results in 

(A-jB)e j at + (A+jB)e“^ a,t ^2A (5.20) 

where A >> B which represents a step function of magnitude 2A. 

A simple example will illustrate these points. In Fig. 5.3 is the 
block diagram of a plant whose transfer function is known. It is required 
to determine the output c(t) for a ramp input of unit slope. 


2 


(3 + 1 )(,s +2) 



C (s) 


Thus 


C(*s) = ■ 


Fig. 5.3 
2 

s^(s+l) (s+2) 


’£5.21) 


By partial fraction expansion and the inverse Laplace transform, one can 
obtain the exact solution 

c(t) = t - | + 2e -t - A e~ 2t (5.22) 


In Fig. 5.4 is the equivalent circuit model for the plant given in Fig. 
5.3 for computer analysis with NASAB. 



Fig. 5.4: Circuit Model for Fig. 5.3 

VI = input r(t) V2 = VC1 ;_ V3 = -VC2"= output c(t) 

Desired Transfer Function W3/W1. 
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Since the raanp response is desired, the HASAP program will evaluate the 
residues of the poles of the function 


C(s) = (i+Jei)(s-te)(»+l)(«+2) 


Res(s 

--1) = 

:^pW2 



l+a 

Res(s 

= -2) = 

-2 1 
O ~ ”* o 

W 2 2 

Res(s 

= d a) 

-d - - d 

" Pf(l+J«)(2+J«) 2 et 

Res(s 

= -da) 

3 i 

" ®(l+3«)(2+d«) “ 2 a 


The computer results for the transient analysis of the circuit in Fig. 
p.k are given in Fig. 5*!?* Rote that for # = 0.001, the coefficient A and B 
are -0.75 and 500 respectively. Indeed B > > A and thus by use of (5.18) we 
obtain 

c(t) a* t-1.5 for large t 

which agrees with the -exact results obtained from (5.22). 


• NASA P RAMP RESPONSE 



NONE 
RAMP 

RESPONSE 


VI 1 
R1 2 

2 1 
3 1 

Cl 3 
V 2 1 

1 IF 
4 1 VC1 

RE 4 

5 0.5 

C2 5 

1 IF 


V3 1 6 1 VC2 

OUTPUT 

VV3/VV1 

nm i,5 

EXECUTE 
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TRANSFER FUHCnOW VV3/W1 


HAW 


RESPUW. 


oaoooe-oi 
_ftU 2006-00- 
O*1?0O£ 00 

_ft.jsoqs_M_ 
0+8L00E 00 


O.572Z0M96-O5 

0.6&ri419^6E-04 

0.22O9030E-03 

_JU5Aiiiiyi^0J_ 

0. J 0040990 E-0 2 
_ 0.169ft».iag.E^0.g_ 
0 ,26ji20355Er,02 
~ 10097E-02 








Although a few examples of HA.SAP printout of step response appeared earlier 
in this- jaanual we include a specific example here * This step response is for 
the Fiseriberg control problem shown previously in Fig. 4.3. The NASA? printout 
is shown in Fig. 5*6. We shall refer to this problem again in Chapter VII. 

We are not including any examples of impulse response since all the pre- 
vious references to MSAP only show such printouts. 
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RAMP RESPONSE 



We continue with the Newton Gould and Kaiser problem, see Fig. 4.9, to 
illustrate error responses. In Fig. 5-8 gives the step error. Figure 5.9 gives 
the step response. The latter is included here to emphasize that only ja single 
NASAP computer card need be changed to get the alternative response output. 
Finally in Fig. 5 *10 we show the ramp- error response for this same control 
system. 
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FIGURES OF MERIT BASED ON ERROR SIGML 

The conventional criteria on this basis are the transient performance and the 
steady state performance. These were briefly summarized in Chapter IV. 

The overall response of a control system is determined by the poles and zeros 
of its transfer function G f (s) - (s) / D f (s); it is not possible to tell it 

from the degree of (s) and the difference of the degrees of D f (s) and N f (s) 
alone. However it is found that for the same an optimal system with a 
smaller SD f is better than such a system with a larger SD f ; where we have let 
5( ) denote the degree of the polynomials. Similarly for the same &D f , an op- 
timal system with a smaller ( 5D^ - 6N f ) is better than such a system with a 
larger (6D^ - ) . Therefore, given a plant, for which we can choose the degree 

of 6 d^ and §N f , it is desirable to choose a smaller and a smaller ( SD^ - 6N^) . 
How the smallest possible ( SD f - 6tT f ) is. governed by the given plant. 

It can be shown [CH l] that the absolute minimum of is 6D - 6N; where 
D and N refer to the uncompensated system. To achieve this minimum, is 
required to be zero. In general > SET, unless some pole-zero cancellations 
are employed in the design. 

The transient performance of, a system may be specified by percentage overshoot, 
rise time and settling time. These specifications are dictated by the poles 
and zeros of the transfer function. Since we do not have control over the 
zeros of a /system, usually we just try to put the poles of the over-all system 
in some desired location. For a second order transfer function with a constant 
numberator, the desired pole locations can be readily determined from the trans- 
ient specifications. For high order transfer functions, the concept of dominant 
poles can often be used. Then a pair of complex conjugate poles is located as in 
the second order transfer function and the rest of the poles are located in the 
far left half plane with real parts at least ten times as large as the real 
parts of the conjugate poles. In choosing these poles, the steady state per- 
formance should be kept in mind. 
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The steady- state performance of a control system G f (s) depends only on the 
coefficients of G f (s). Let 


G.(s) 


h Q + b^s + b^s + 


+ h _ s 
n-1 


n-1 


a Q + a^s + a g s + 


, „ n-1 , n 

+ Vl S + 3 


If the steady- state error due to a step input and 'a ramp input are required to 
be smaller than k a /^, then we need, respectively. 


| - | fi k/100; and a = b , | -1- i | s k/lOO. 

o o 

Hence the steady- state performance of a system can be rather easily controlled, 

When using these criteria it is necessary to check the response of the chosen 

overall transfer function with an analog or a digital computer to be sure that 

it is satisfactory before continuing the design. 

In addition to the conventional criteria just discussed, we shall mention 

a few other criteria based on the error signal that serve to make the control 

system "optimum" in some sense. 

(i) I1AS criterion (integral of time -multiplied absolute -value of error): 
This criterion was first introduced by Graham and Lathrop [GA 1]* For 
a given plant, the problem is to design an overaall -system which mini- 
mizes 

J t|r(t) - c(t) |dt 

o 

where r is the reference or desired signal and c is the output of the 
over-all system. It is clear that |r(t) - c(t) | is the error between 
the desired signal and the actual output. The multiplication of t on 
|r(t) - c(t)| provides an increasingly heavy penalty for a sustained 
error. Using step functions as reference inputs, Graham and Lathrop 
obtained, by analog computer simulation a set of optimal transfer 
functions [GA l] . 
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To choose an optimal transfer function G^(s) for a given- plant 
G(s) note that if all the forward paths from r to c pass through the 
plant, then the zeros of the plant (the roots of N(s)) will he indepen- 
dent of how the compensators are introduced. Consequently if the 
numerator of G^Cs) does not contain all the zeros of G(s), the missing 
zero must he cancelled hy a pole. The advantage of this criterion is 
that it is very selective; however it cannot he studied analytically. 
This criterion is not widely used, because a complete list of optimal 
transfer functions is not available. 

(ii) Quadratic criterion : The optimal system is the one which minimizes 

to 

+ (*(*) - =(t)) 2 3 at 

O ' 

O 

where u is the input to the plant. In this criterion, if u is not 
included, the optimal transfer function will always he unity and the 
required compensators may not he physically realizable; furthermore 
the magnitude of u may he large and the .system will he saturated. 

The optimal transfer function for a given plant and a given reference 
input r can he obtained hy applying Chang's root- square -locus method 
[CH l] as well as hy using the dynamical equation description [see 
AT 1] 

The design hy using the quadratic criterion can he solved ri- 
gorously. However there are three arguments against using this 
criterion. First, it is not very selective [GR l]. Second, the 
criterion is chosen mainly for mathematical convenience rather than 
practical reasons. Finally and most seriously, the resulting optimal 
transfer function may not he realizable in practice. If all the 
zeros of the transfer function of the plant G(s) = have negative 

real parts and if the reference input is a step function, the optimal 
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transfer function is of the form and 6D f (s) = 6D. In this 

case, it is easily realized. However for plants with nonminimum- 
phase transfer function, the design by using quadratic criterion 
might have difficulty. Another difficulty in using this criterion 
occurs when the reference input is a ramp. 

(iii) ISE criterion ' (the integral of the error squared): Here one seeks 

to minimize 

J e (t) dt. 
o 

To obtain any of these criteria from UASAP use of an external integration 
subroutine is required. To illustrate this we follow Beck [BE 1] and 
use the integral of the squared error as the performance index although 
other criteria can be applied with equal ease since an analytical so- 
lution is not required. The formation of the chosen performance index 
is indicated in Pig. 5.^ taken from [BE l] where it is shown that the 


Input 



To 6^) Converter & 
Digital Computer 


Fig. 5*?: Index of performance formation 

ISE is generated in the analog computer and returned to the digital 
computer through the analog to digital converters. In the hybrid 
computer application an optimization algorithm operates upon this 
output . 
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CHAPTER VI 


SENSITIVITY ANALYSIS 
VIA INTRODUCTION TO SENSITIVITY 

The NASAP sensitivity results, given in tabular and graphical form, 

-can be used to predict the percent change, absolute change, and modified 
value of the transfer function to changes in a particular network parameter. 
By definition the sensitivity of some real function* o (in NASAP, ReH, timH 
| H j and $) to a change in a real parameter x (in NASAP, resistance, capaci- 
tance, inductance and dependency value) is defined as 


C ° A x *« 

s “'x ” a dx 


d a 

_a _ dlnva 

dx dlnx 
x 


( 6 . 1 ) 


By rearranging (6.1), the differential da can be related to the differential 


dx. 


d a=‘ 



dx 

x 


( 6 . 2 ) 


For an incremental change. Ax, in the parameter x the incremental change, 
A a , in a can be approximated from (6.2) 

^ a 


A* dc 


0 S a i 


Ax 

X X- 


(6.3 ^l) 


If the change in x is expressed as some percentage of x, then (6.3a) can 


be expressed as 


A at 


C a l£_ 

a O x inn 


(6.3b) 


where Px is the percent change in the parameter x. The percent change in- a 
is easily found by dividing both sides of (6.3&) by‘ ot , thus 
a 

Pa^5 x ^ (6.4) 
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where P .e*. is the percent change in the function a . The modified value of 
' v , , ^called* a T , can be expressed as 



a 1 A i a+ Ace 

(6.5) 


■ o'*a(l +S X 

(6.6a) 

or 

a ' o(i + S x fgs) 

(6.6b) 


Naturally the smaller A x or the more accurate- 9 ' a 1 s andiP<£i' 
will be. The accuracy will also be enhanced if a is roughly a linear - 
function of x over the range of x under investigation. 

If indeed’ a * is a linear function of x of the form 

• a= Kx, (6.7) 

then from (6.1) it is seen that 

( 6 . 8 ) 

^x *' 

Furthermore from (6.1) if a is independent of x, then 

C “ T. 0 (6.9) 

'— > x 

or the function a is insensitive to changes in x. 

Similarly the root sensitivity printed out by NASAP can be used to 
predict the new pole and zero locations for a small change in some network 
parameter. The root sensitivity is defined as 

. 15 1 



where p is a zero or pole of the given transfer function. The differential 
dp can be expressed in terms of the root sensitivity as 


dp 



dx 

x 


o P' dlnx 
x 


C6.ll) 


Note that dp will be complex since tr? is also a complex number while 
dx . 

— is a real quantity. 

The incremental change in p, Ap , f 0 f' an incremental change in x is 
derived from (6.11) i 


Ap 



( 6 . 12 ) 


If Ax is expressed as some percentage of x then (6.12) becomes 


Ap 


oV 

x 100 


(6.13) 


where Px>is the percent change in x. 

VTB .DERIVATION 0E SENSITIVITY FORMULAS 

Although the formulas used in NASAP to calculate the transfer function 
sensitivity were derived from a tagging technique on the loops of the flow- 
graph (see l]),j these sensitivity formulas can be obtained by using 
simple calculus. Suppose the transfer function is given 


H(s) 


N(s) 

D(s) 


(6.14) 


where N(s) and D(s) are polynomials in s. 

Since H(s) is the transfer function of a linear circuit it can be ex- 
pressed as a bilinear function of any element in the circuit. That is, both 



N(s) and D(s) can be expressed as the sum of two polynomials in s where 
one polynomial does and the other does not contain the specified element. 
Thus 


H(s) = 


N(s) 

^(s) 


A(s) + x B(s) 
C(s) + x D(s) 


(6.15) 


where x is the specified element in the circuit and A,B,C, and D are 
polynomials in s. 

The sensitivity of the transfer function H(s) to some parameter x 
is defined as 


$ H(S) A dH/H = dlnH 
x =^dx/x dlnx 


(6.16) 


Thus to find the sensitivity of H(s) to changes in x in (6.15), one 
obtains 






(6.17) 


but 


n n(s) 

O _ x d N(s) 

x “ N(s) dx 

- ss -t < A+Bx) 

_ Bx 
A+Bx 

Similarly 


S° iW 


Dx 

C+Dx 


Substituting (6.18) and (6.19) into (6.17) yields 
H(s) 


5 


Bx 

A+Bx 


Dx 

C+Dx 


(6.18) 


(6.19) 


( 6 . 20 ) 
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This expression remains unchanged if 1 is added and subtracted from 
the righthand side. Then we can write 



_ Bx A+Bx _ Dx C+Dx 

A+Bx " A+Bx C+Dx C+Dx 

A_ + _ Q — 

A+Bx C+Dx 


A(s) C(s) (6,. 21) 

N(s) + D 1 (s) 


This equation can also be derived from the tagging technique on the loops 
of the flowgraph and is the one used in subroutine SENSS of NASAP to deter- 
mine the transfer function sensitivity. 

VIC "discussion OF SENS ITIVTTY ..formic sens- 

In subroutine SENS of the NASAP program* are calculated the sensitivity 

s, 


expressions 


k ReH, 


ImH, 

O* 



where H(jm) = ReH(j(o)+jIm H(jai) 
and H(jto) = |H(jw)|e^ . 

Again the tagging techniques of 1] are used to determine the sensitivity 
expressions. The basis for the tagging procedure is that the transfer 
function H(s) can be written as in (6.15) where x is the sensitivity para- 
meter.. Since the polynomials A, B, C, and D are, in general, complex quanti- 
ties for s = jw, (6.15) can be rewritten as 


\ _ (Re A + x Re B) + j (Im A + x Im B) 
^ (Re C + x Re D) + j (Im C + x Im D) 


( 6 . 22 ) 


where A(jo>) = ReA(jw) f j Im A Cjw) etc. 

After some mathematical manipulation, the right side of (6.22) can be 
separated into its real i |md imaginary parts..., s Thus_ - ^ 



o' , ?R (ReA ReC + ImA ImC) + x(ReB ReC + ReA ReD + ImB ImC + ImA ImD) + x^(ReB ReD + ImB ImD) 

ReH(ja>>. -i— = « 9 7 2 (6.23) 

U [(ReC) + (ImC) + 2x(ReC ReD + ImC ImD) + xH (ReD)*+ (ImD)* 


and 
ImH(jio) 


_ N I _ (ImA ReC - ReA ImC) + x(ImA ReD - ReA ImD + ImB ReC - ReB:. ImC) + x 2 (ImB ReD - ReB ImD) ^ ^ 


E< 


(ReC) 2 + (ImC) 2 + 2x(ReC ReD + ImC ImD) + x 2 (ReD) 2 + (ImD) 2 ] 


Thus by use of (6.17) one obtains 


S 


ReH 


■S R - 5 


= *_ j. n .ii B 

dx n R Ddx" 
K 


ReH 


dD 


S._ N^D dx N R N R dx 


(6.25) 


dML 

wftSre r= — = (ReB ReC + ReA ReD + ImB ImC +• ImA ImD) + 2x(ReB ReD + ImB ImD) 
dx 

and = 2 (ReC ReD + ImC ImD) + 2x j^ReD) 2 + 


(6.26) 

(6.27) 


and N_ and D are defined in (6.23), 
R 


By a similar use of (6.17) one arrives at the expression 



x r* i ■ dD 

N D \dx ~ I dx. 

I ' J 


where — ^ = (ImA ReD - ReA ImD + ImB ReC - ReB ImC) + 2x(ImB ReD - ReB ImD) (6.25 

and is given, (6.27), and and D are defined in (6.24). 

Equations (6.25) and (6.28) are used in subroutine SENS to evaluate the 
sensitivities of the real part and imaginary part of the transfer function to 


changes in the parameter x. 


The sensitivities 


ies \ and 


are evaluated in terms of the sensiti- 


vities found in (6.25) and (6.28). By definition 


|h| = V( ReH > 2 + Cm) 2 


Differentiating this with respect to the sensitivity parameter x yields 

_ „ dReH , T „ dlmH 
ReH — - — + ImH — — 

, dx ' * dx 

4-1 H| = (6 

\ dx 1 1 

V (ReH) ^ + (ImH 


By definition 


o 1 1 X 

0 X = w 


Substituting (6.30) and (6.31) into (6.32) results in 

I „ | „ „ dReH , _ dlmH 

O 1 dx dx 


Recalling that 


0 reH A _x_ dReH 
* ReH dx 


^ImH ^ x dlmH 
V '^x ImH dx 5 


one can simplify (6.33) to 



(6.36) 


q|h| (ReH) 2 $ ReH + (ImH) 2 

Me = X -ThI 2 


The sensitivity of the phase of the transfer function to changes in x, 
is also easily obtained. By definition 


By implicit differentiation of (6.37) with respect to x, one obtains 
Tj dlmH • _ „ dReH 

2. d* ReH — " M IT 


A _x d£_ 
x dx 

one obtains by substituting (6.38) into (6.39) 


. B _ dlmH _ dReH. 

-rr xReH — x-XmH— 

v 2 l dx * i ». dx , - 


By use of the definitions in (6.34) and (6.35) and the relation 


ReH = H cos 


the expression in (6.40) can be simplified to 


■vlmH ^ReH 


S * = ;iM ReHM )( S * "S 


Equations (6.36) and (6.42) are used in subroutine SENS to determine 


^ ^and , respectively. 

This is possible since the quantities <f> , |h| , ReH, and ImH have been 
previously calculated and stored during generation of the Bode tables and 
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ReH ImH 

plots while O x and ^x ^ ave been determined earlier in subroutine 
SENS . 

V3D DISCUSSION OF REVISED SENSITIVITY SimOUTINE 

In Appendix B is a revised version of subroutine SENSS which does the 

calculations of SENS and SENSS in a more efficient manner. This version of 

SENSS requires only 3/4 of the core storage required by the present SENS and 

SENSS and it utilizes a simpler algorithm that greatly reduces the number 

of mathematical operations required. This saves execution time and should 

increase the accuracy of the sensitivity calculations. 

The version of SENSS given in Appendix B uses the same tagging procedure 

that is used in the present SENSS and calculates with the use of 

equation (6.25) as does the present SENSS. However, the outputs of the 
c [H| cf> c ReH Q ImH 

present SENS, that is jO x j andw^ are related to the 

real and imaginary parts ofg H ^ W ^ Thus the rather complication sensitivity 

x 

expressions now used in the present SENS and described in equations (6.23) 
through (6.42) are completely avoided. 

c H(jw) 

Since, in general, the sensitivity expression O is a complex 

quantity, it can be written as 


s H Uw) „ ReS h( 3*) + jIm s 

'x X X 

But 

H(jto) - | H(joj) 1 e 3 ’’ 

where 4 is the phase of the transfer function H(joj). 

Thus one has . , 

HUco) |H|e^ 

O X "Ox 


'H(jm) 


=s 11 + s 

■V ^ 




(6.43) 

(6.44) 


(6.45) 


1 64 



Let us examine the rightmost term of (6.45) more closely. By definition 


s. 


e*^ A :-,x d j<J> 
x " e# dx 6 ‘ 


(6.46) 


By use of the differentiation chain rule, the derivative expression in 
(6.46) can be simplified to 

e j*f> = ^ M. 


dx 


dx 


Substituting (6.46) and (6.47) into (6.45) yields 

H| 


dx 


(6-;47) 


(6.48) 


Using the definition given in (6.39) this equation can be rewritten as 


S 




\ 


s* 


(6.49) 


If the sensitivity parameter is a real quantity, then the expressions 

and ^ ^ ,,will also be real. Thus if (6.49) is compared with (6.43) 
and the real and imaginary parts equated (under the assumption that x is real) 


then one obtains 


and 


S H = Re S 


HO) 


>S - ii»S 


H(joi) 


(6.50) 


(6.51a) 


s; 

where x is a real variable. 


H(jw) 


(6.51b) 
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ReH 


The sensitivity of the real part of H(SfV) to changes in x, 5 can' 

X -»H(ju) a 

'jc 

By definition 1 


also be obtained in terms of the real and imaginary parts of 

ReH ( 3 co) = | H(jco) | cos ’ (6.52) 

g ReH ^|h|cos <f> 

X X 

igfH| '£ cos $ 


Thus 


(6.53) 


But 


s 


cos^cj?, dx 
deft 

-X tan $ ^ 


(6.54) 


Equation (6.54) can be rewritten as 

S COS (f» Q <f> 

= -({> tan (ft O 
x x 

by use of the definition in (6.39). Thus equation (6.53) becomes* after 
the substitution of (6.55) ,> 

ReH 


(6.55) 


S - s 


H 


tan ■ 


s \ 


(6.56) 


However, after (6.50), (6.51a) and (6.37) are substituted into (6.56), 


the expression becomes 
ReH 


S - HeS. 


H(jco) 


ReH 


Im- 


$ 


H(jto) 


(6.57) 


A similar expression can be derived for the sensitivity of the imaginary 

, ImH 

part of H(jw) with respect to changes in x, b . By definition 


ImH(j(o) = |H(jco)| sin (ft 


(6.58) 
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Thus 


ImH 


H sin ' 


|H| g sin $ . 
x + x 

By a similar mathematical technique, it can be shown that 
sin <f> , _ 4 


S, 


tan i 


Substituting this expression into (6.59) gives 

■ S * 


5 ImIi _ 5 i H i 4, 

x x + tan i 


(6.59) 


(6.60) 


(6.61) 


which can be further simplified by the substitution of (6.50), (6.51a) and 
(6.37) 


S “ = ReS Im S H(j “ } 


(6.62) 


Equations (6.50), (6.51b), (6.57) and (6.62) are used in the version of 
SENSS given in Appendix B. This is possible since the real quantities , 

ReH, and ImH have been calculated earlier in subroutine BODE while the complex 

HU“) 

quantity S x is calculated in the Appendix B version of SENSS. 

Note th^t (6.57) involves only 3 arithmetic operations while (6.25) 
involves 8 arithmetic operations plus the numerous operations involved in 
equations (6.23) and (6.26). The same comparison can be made between (6v62) 

. c | h| 

and (6.28). With regard.' ;to O , (6.50) involves no arithmetic operations 

- c'+ 

while (6.36) uses 8 operations. Similarly in determining equation 

(6.51b) requires one arithmetic operation while (6.42) involves 6 operations. 
Since the sensitivity calculations must be redone for each frequency value, 
the number of arithmetic operations is quite substantially reduced with the 
version of SENSS given in Appendix B. 
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VIE ROOT SENSITIVITY 

The sensitivites of the poles and zeros of the transfer function 
to changes in some specified network parameter are also determined by NASAP 
by use of the tagging technique and the formulas given in [KU1, PAl] which 
will be summarized here. 

Suppose we have a polynomial in the complex frequency variable s, 

P(s), which can also be expressed as 

N 

P(s) = 0C( s) + x{5(s) = £ a, s^ (6.63) 

k— o * 

V' 

where the degree of P(s) is N and P(s) has only positive powers of s and 
x is the sensitivity parameter. Furthermore the roots of the above polynomial 
are known, i.e., N values of s are known such that 


P(rJ =0 i = 1,2, • * * N (6.64) 

where r^ is the i^ root of P(s) . 

ri 

If the order of the root r^ is m, then the root sensitivity, s^ , as 

defined by , 

ri dr i 

s r Ax -z — (6.65) 

x dx 

can be expressed in terms of the given polynomials- ct fcs-J and $(s) 


,m-l 


ds 

d^_ 

ds 1 " 


m-1 


(x3(s)) 


(s) + x3(s)) 


If r. is a simple root of P(s) , then 


i~l,2, * * * N (6.66 y- 


T. 

i 


x$(r J 


ds 


(oc(s) + x3(s)) 


(6.67) 
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where M is the number of distinct zeros and A(s) , B(s) , and N(s) are 
defined in (6.23). 

Similarly, if is a simple pole of H(s) then the pole sensitivity with 
respect to x can be written as 

V ± xD(P ± ) -xV(-p ± ) 

X ~ (C(s) + xD(s)) I “ I»i(s) 

S=P i 

where N is the number of distinct poles of H(s) and C(s), D(s), and D^(s) are 
defined in (6.68). 

Since the polynomials A(s) , xB(s), C(s), and xD(s) have been determined 
by the tagging process during the evaluation of the loops in subroutines 
FLGRPH and HIGORL, the sensitivities of the poles and zeros are easily obtained 
by differentiating the denominator and numberator polynomials respectively 
and by evaluating the resultant polynomial and the appropriate tagged polynomial 
at the given pole or zero. These calculations are performed in subroutine 
ROOTSS. 
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VIF EXA 1 MHES 

In the foregoing sections of this chapter we have indicated many possi- 
bilities for sensitivity analysis with the aid of HASAP. We shall illustrate 
a few of these as follows: 

1 * Taking advantage of the background developed in Chapters r V and V 
on the unity feedback control system with lead cascade compensation see Fig. 

6.1a (also Figs . 3 .16 and 4 . 7 ) we obtain the sensitivity of the transfer func- 
tion W 4 /W 5 with respect to resistor Rl. This is a judicious choice of circuit 
element since one of the time constants of the control system plant is 

\ - (El) (Cl) 

Hence if Cl is constant in value the sensitivity determined for Rl is the same 
as that for the time constant 1^. The WASAP printout for this example is shorn 
in Fig. 6 . 1 . The corresponding zero and pole sensitivities are given in Fig. 

6 . 2 . 

2 . The second example Fig. 6.2 follows up with the uncompensated unity feed- 
back control system for which we obtain the ramp response in Fig. 6.3 and the sen- 

H 

sitivity function . Thus we seek the sensitivity of the transfer function 
W 4 /W 5 with respect to system gain K, in this case K = O.83. The KASAP printout 
is shown in Fig. 6.4 with the corresponding zero and pole sensitivities given in 
Fig. 6.5. 

3 . The third example furnishes sensitivity data, specifically S H for the 

K 

multiloop feedback compensated control system with K = 38. Since the pertinent 
HASAP model and transfer function print out were given in Fig. 5 , they are not 
repeated here. The corresponding gain sensitivity print outs are given in Fig. 

6.6 and the zero-pole sensitivity in Fig. 6.7. 
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4 . For the final example we use one of the control subsystem of the 
Ranger space vehicle taken from Dorf [D 01 , Problem 4 . 5 ]- The mission in 1965 
was to scan the lunar surface with TV and other sensors. The requirement on 
the altitude control subsystem was to stabilize and control the Ranger space- 
craft from second stage separation to lunar Impact. Briefly a high gain an- 
tenna furnished input signals to the earth horizon sensor which in turn fed the 
gyro loop. The latter consisted of the spacecraft, as the plant, with the 
altitude gyro in cascade as shown in Fig. 6 . 8 z along with its HA.SAP model. To 
give some feel for the cut and try process in determining the gyro loop gain 

K that will keep the step response overshoot under the required 5 $, we show 

H P-4 

for three values of gain the print outs of S^, zero-pole sensitivities S„ 

K H 

and the corresponding step response. For convenience in examining the print 
outs we tabulate the pertinent figure numbers. 


K 

step 

response 

s H 

s k 

s P i 

S H 

6 

/ 6.8 

Fig. 6.9 

6.10 

20 

6.11 

Fig. 6.12 

6.13 

75 

6.14 

Fig. 6.15 

6.16 
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0. 60019) lb 00 
. 11. f. 7 I ,iV.2E (.!■■ 
O. 725HM'‘»b 00 

LL«_6l/tJ£ 0 -lU 1. _0 0 — 

U. 0090001b 00 


0.3/56453t 00 
_U,655i2135b^l- 

-0.3H02JO7E 00 
-U. 1307276b 01 


-0. 2734890b ('1 
-» .43vxrt».yiT iiv * 
-U.473U998I:-<’1 


-0 . 10666 j 9fc 01 
-fi.inv» ; /3h (11.. 
-0. 1 0'*0956b 01 


.>1,1.6 ( iii. 6 I.IH 

0. 12030906 00 
0.1676536b .00 „ 
0. l66|.H6E 00 
fr ■ t H? mir .10— 

v. 1 95 1 ;oxi; o > 

_JU. 201.22311; _UU._ 
0,1939622b 00 


-0.81229156-01 
rO.256Mjy2fc.00_ 
-0.4734553b 00 
- 0. . Jm 411.4 0 L_0 0— 
-0. 5379563b 00 
i lfifiV4 2t- 01 
-O.j 2039/91: 01 
-O. I239»~49g 01 
-0 , J 234407b 01 

IA ffl LUL 
-0. 1172447b 01 
~0..137?2Hfe 01- 
-0. 1 105636b 01 
-0 . I 079099b 01 
-0.1057704b 01 
■riL,J.O'*L2AU_0.1_ 
-0. 102876 7b OJ 


-0.1005313b 01 


M .450 .11 ct) J 
0.33679626 00 
„ 0.3V52041L 00 
0.45624616 00 
n.. i sz/.A3i7iF_ nn _ 
0.592) 40?r 00 

0 .650 7 504 1 OU. 

0.7215140E 00 

p . /7 a9 B6ZE_C 0— 

0.02085356 00 

n.<l47'iSI7f n»__ 

0.851 11096 00 

0, 82o9603E .00— 

0.75405246 00 
__0. 6538234E 00„ 
0, 53063378 00 

n./.04i^H4g -in 

0.29480956 00 
_ CL. . 2.0. L (16 9th._U.0_ 
0.14469886 00 

OjJJI lU3i Zi — CIQ_ 

0.71379136-01 
0. SIfii'2 5 5r-(H 
-0.4 1467 101. -01 
-0.34448956-01.- 
-0.295776^-01 
_^26 i9349Et'J 1 - 
-0.23000976-01 
-o. 72o-)q48r-;n__ 
-0.20655416-0) 
_rD..J.2:UBi9L=01_ 
-0 , 10400596-0 l 
_^Iu237jQ'32Erai_ 
-ti.1632749r.-01 


-U.9860165E 00 


-O.i 00 1859b OL -t,.1632749r-01 
-(i . , n 00994 h OJ -O. I 5t 


-0.1000466b 01 -0. 14?571)6-01 « 

-|> . "■ 0001 64. C._03 “0...1323J.0O ts01_ 

-0,9999768c 00 -0 . L222457t-01 

-0. 9998509b UO -0 . 103 11 74fc-01 

-P.699fe4V0b UO -n.94yl5 7.TF-OX... 

-O.V998489b 00 -0.05626446-02 



.UO ( it: Si U li LH JJ jL_|,t!bjL5 L ' l5.( ALti i il LU_LU6£iUJ6 ( ,**F 1 (11).) >. 
-0.34269446 00 -0,«V146926 00 -0.472o32ftE 00 


-0.72171046 00 

sl-’.l UUx\lY ,)» 




-0.//967616 00 


-0.33690126 00 

— ~i>« 2i?i , ?iV7F nn 


-y. t jjAyt rt c»i 
=u.aQalHi9fc on 


ikV«j‘i/u yy 

uq 


-U .///4744r O o 
-n.u n>79ip nn 


•0.9 10^2116-01 


-U.1336/Q26 00 
1/ t?(,F f t? 


-0.71232 75E‘ UO 
-b.767<s5; te t»n 


-0.9ybg?40t 00 

-ft. i -iVfi U?h 


-0.1417532? 00 
■-0. 109.Stt64E-flO- 


-0.8573431C-01 
-■iU71H3lV«lf H ni 


-0.42522206 00 -Q.l090?87fc 01 

■aU. HUHOIF (II -Oai-9l95-/. (l E_UU-_ 

-0.41V94256 CO -0.324 7210fc 00 

-Ji0.6b!]7045lU301 rO. l46056lt 00. 

0.ll03o$9E 00 -0.2/817376-01 

uUSZyiH)- .cn 


“0.7052439r-01 

UEdU_ 


-0.1221380C 00 
t 0.1045395E 00 - 
-0.2/520S1E 00 
=&-lS,739ftAF flQ 


• 273<Jb-66 00 O.C»6i9046-OJ 


0.44010416 CO 
JA42.74 7h . 


0.o5043236-6l 
-a« 55847 •> it? -ui 


-0.53C45B5G 00 
={U66JL7912E-0U_ 
-0.5389352E OO 
L..95.554-40t_00_ 


“0.1146423E 0) 

J?ii7^7»P r. . 


-0.13243m Cl 
LLiiijAiL. f.n 


.43012046-01 


-0.8457/766-01 U . fc43V702t-01 

-U . 3403790b -0 l 0 . 1 23 1687L-01 

t P.'.4*7?>4)-’»-/ 


-0. l3l)229*)E Ul 
_r0..JL462824E_0JL_ 
-0.13290J6E 01 
_=Q..15i;i805£..0JL_ 
-0.(6734056 01 
_r(L«.l 6ti&79 tl E_fll. 


•0.2u006/ Jb-y 1 

•U. 


•0.13351246-01 

LQ J JJ31^365t.-01 a, 

-U.VbV2j76t--(<2 0, 

Cb 


.•»y/ //346-02 
1 /««-- o/ 


23012026-02 
a3Ski23±-u<L_ 
. J0647»*66-03 
432l*40 e -tn 


-0.1684965? 01 
-=lL.±I10.446t_OJL_ 


-0. 173516*E 01 
_rO. 1.760.4 4 2t_01_ 
-0.(7b7o8or ni 
-fl.mii44^P >n 


-0.61168356-02 


-0.40422606-02 -0 

■ Ol 47f)/459 k-C/ -0 


>0*2632 77.26-02 


2024&46e-03 

,91396136-03 

.-4-6306906-04 


64 1 1 1 516-04 
:»H60373b-»4 


-0.1«45967F 01 
_rO. 18784076 01... 
-0.1VU766F 01 
-=a.a2AA9oaf..oi_ 


-0.1V36667E 01 

-f). 50750 ZAP Ml 


-U.1696406E-02 -0.65625866-04 -0.2O66379E 01 


e 

c 

0 

8 

0 

O 

© 


0 

“ © 

— O 
“ © 

.. © 


.© 

© 

© 


o 


_© 

© 














■6LX 


© 


NASAP-PRUBLEil -UNCOMPENSATED. PLANT: # 




He*rt 

• 

-RAHP 

RES-ON.Se 


VI 

1 

2 

C.B3 . R3 


LI 

2 

1 

l.H 


V2 

1 

3 

1.0 111 

1— 

Rl 

3 

4 

1. - - - 


Cl 

A 

1 

1 r 


V3 

1 

5 

1.0 Vvl - 


R2 

5 

6 

0,2 • 


C2 

6 

1 

X f - 


V4 

1 

7 

1.0 Vv.2 

J 

-R3 

6 

7 

-1-.0 


V5 

1 

B 

1.0 


OUTPUT 



VV4/VV3/’. 1 
FREQ 0.0 1.0 0.05 
TIME Si>. 

EXECUTE 



© 

© 


© 

© 


© 

TRANSFER FUNCTION W4/VV 5/VI 

" ” *•-* 

' NUMBER Or LOOPS PER URDER ■ 
1= 4 

© 

© 

{ ) 



~ © 

© 

i H(S)*- *t. 150c 00* 



.... © 

© 

{ 4.15E 00 +5.00E 00 S 

+6.00E 00 S +1.00E 00 S ) 


© 

o 

„ 7 con nc TR P’S , PnMfTTfiN . » . - >» » — 




© 


i * “ ^ r 4 uKU Ur | [\OnjrCA ruiw • X lhi ~ 




o 




© 

o 

a 

POLE OF TRANSFER FUNCTION 
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CHAPTER VII 


SPECIAL CONTROL SYSTEM EXAMPLES 

In the previous chapters of this manual we have used a variety of control 
circuit problems to illustrate various facets of computer aided design with 
NASAP. The examples were basic to control theory with e:xplicit or implicit 
connection to aerospace applications. In this final chapter we have selected 
a few control system examples that emphasize the capabilities and also the 
limitations of the present version of NASAP as modified at the Moore School. 
From these examples it will be apparent that, in general, NASAP can be an 
effective aid to the control system engineer dealing with moderately complex 
linear control systems. 


210 



VI IA EXAMHH I3W0LVING TIME DOMAIN APIftOXIMATION 


Approximation of system response is a long standing problem related to 
control system analysis and synthesis. It is often necessary to approximate 
a system before suitable compensation can be applied. Eisenberg [El 1] has 
described a technique for approximately identifying high order system responses 
based on (.certain characteristic responses to a unit step input. Many feedback 
systems exhibit time domain response to a unit step input that includes an 
overshoot followed by variations which subsequently settle to a steady state 
value. The approximate system response is generated by a closed loop trans- 
fer function whose open loop transfer function contains a "transport lag and 
a first order lag. This technique only requires knowledge of the first peak 
overshoot, the time to peak and the settling time of the unkn own system 
response. 

We represent the unknown response by that of a unity- feedback control 
system with G(s) representing the plant. The general form of G(s) is assumed 
to be 


G(s) 


K e 
P 


*sTd 


TS + 1 


(7.1) 


This particular approximation is quite convenient and yields the system transfer 
function 


G(s) 

K S " 1 + G(s) 


5=i -sM 


i + J2. P - S ® 


( 7 - 2 ) 


To facilitate further consideration of this expression, we normalize it with 
respect to T^. Defining 


s = TqS = crT^ + jcoT, = cr + 
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we can write 


|(s) =■ 


a e 


+ p + «e 


(7-3) 


where 

K T , 

* = p = — * * = V 

The problem of system approximation is now reduced to that of determining 
the parameters and p such that the denominator of (7.3) has a pair of clearly 
dominant complex conjugate roots which will exhibit a second-order response. 

In particular it is desired that the responses shown in Fig. 7.1 he similar , 
for the set of conditions 


M ~ 14 = M, T , ~ T = T, T ~ T = T (7.4) 

ns' ph ps 9 sh ss s w ' 

where the values of T^j and are measured from the unknown system step 
response. 

One way of accomplishing the selection of a and £ is to use the generalized 

curves found, in [El l]. An alternative is to find a suitable electric circuit 

model for G(s) in (7.1) and then use NASAP. The new feature here is the ex- 
— sTd 

ponential, e . We can obtain a rational function representation for this 
exponential by a Fade approximation. Specifically we use the biquadratic 

2 

s - as + b 
2 

s + as + b 

In Chapter III* we showed how to obtain a ladder network for this with 
negative as well as positive elements. The input impe'danee of the network is 
the desired circuit model. 

To- illustrate this approximation technique we again consider the seventh 
order system with an open loop transfer function 
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C b> Second-Order System Response 
Deloyed in Time 


Fig. 7*1 System Approximation by Time Domain Response 
(Subscripts M h tl and "s" refer to higher order 
and second orders system respectively) 


G(s) ^ 

(6s + l) (2s + lA (s + 1 a 


(7-5) 


The step response for this in a unity feedback system was given in Chapter V 
as the "Eiseriberg” problem. Consideration of the actual closed-loop.-, unit 
step response of this system, the lower response curve shorn in Fig. 7.2 from 
[El 1] ; yields the three data point = 31.7$, T^ = 20s, and T gh = 55.4s 2 . 
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APPROXIMATE RESPONSE 


1 

o 

t- 

3 

CL 

t- 

Z> 

o 



V. M- ERROR * I IS % 
T p ERROR = 10 5 V„ 
T s ERROR = 15 V. 


70 


80 90 


• TIME sec. . 


Fig. 7*2 Comparison of System Time Responses. 


To obtain the approximate, system parameters we must fall back on the .well known 
'^relationships for a second-order system 



(7.6) 


T 

ps 



(7.7) 


T 


ss 



(for a 2 % deviation of the response 
envelope from the steady- state value) 


(7.8) 


The determination of £ and o> n and ultimately location of the dominant 
roots involves use of (7*^0 • Since M g , where is a measured quantity, 

we have the numerical value Hence the value of £-can either be computed 

from (7.6) as 



(7.9) 


€ = 


1idm| 

\/ (lnM) 2 + 1 


or it can be obtained from a universal curve of M versus Next to determine 
co . Note from Fig. 7*1 a and b and (7.4) that 


\ = T sh - T ph = T - - T - 


ss ps 


(7.10) 


Substituting the required values and solving for co n gives 


v \> r , 


(7. XI) 


Tg'V 1 - r 

Since £ is known co^ can be expressed in terms of measured quantities as 


= V(XnM) 2 + it 2 


|TnMf 


- l 


( 7 . 12 ) 


Next we find as a function of the measured parameters. From Fig. 7* la and b 


T d T ph “ T ps 


(7.13) 


where T - is a measured quantity and T can be obtained from (7*7). Substi- 
pn ps 

tuting (7*7) into (7.13) gives 


T d = 


- llnMlT sh 

4 - |lnM| 


(7.14) 


so that the values of a> n , and are completely specified from the key 
characteristics of the unknown system response. Thus the values of the 
approximate system parameters are found to be = 5* 7s, t = 89.7s, and 
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Kp = 13 ;o l>. These give an approximate open- loop transfer function of 

( 7 - 15 ) 

The closed-loop response corresponding to this G(s) obtained in the re- 
ference is shown in the Tipper response curve in Fig. 7.2. The percentage 
differences between the indicated characteristics of the two responses are 
also shown. It should be pointed out that altough the responses match quite 

closely with respect to M, T • and T the steady- state levels will not match. 

p s 

This mismatch is unimportant because the steady- state levels can always be 
matched exactly by simply adding gain outside the feedback loop. Also the 
gain of the unknown system response is easily computed from the me asurable 
steady-state level, (c gs (t)), via the relationship 

% s (t) 

The NASAP print out and step response for (7.15) is shown in Figs.. 7.3 and 

7 . 4 . 

It should be noted for aerospace control applications a distinct advantage 
of the technique is that it is not necessary to open any feedback loops to 
effect the approximation. Thus, the system response can be approximated 
without interrupting normal operation. This is in contrast to some methods 
of system approximation, this technique is not hampered by the rare cases 
where the open-loop response is oscillatory. There are other methods avail- 
able for system approximation that can be used on a closed-loop basis but 
these methods usually require that the system gain be increased until critical 
cycling is attained (oscillations). For many practical applications, how- 
ever, it is not advisable to bring a control system to the verge of instability. 
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VIIB EXAMPLE OF A CONTROL SYSTEM WITH 'TRANSPORT LAG 

As an ^rtensionj'of the concepts introduced in the previous example, we con- 
sider a given control system that has -a plant with transport lag along with 
three simple lags shown in Fig. 7.5. Lupfer and Oglesby [LU 1] discussed such 
a problem wherein the object was to find a proportional- integral controller. 

They treated this problem with the aid of an analog computer. An alternative 
solution was presented by Eisenberg [El 2] using the parameter plane approach 
to develop a graphical .technique . 



Fig. 7.5 Control system with transport lag 
Reactor process from [LU 1] 


To use KASAP we must obtain a circuit model for the control system with trans- 
port lag. From the previous approximation problem we know how to handle e” sT;D 

by using a Pade approximation that yields a rational function trealized^-by'a' ladder ' 

"r — “ 

network. For the controller of this system we require a different type of 
model. We can write 


K c (T r S + 1) 


( x + rsr) 


( 7 : 17 ) , 


use the model shown in Fig. 7*6 where V7 and 12 are dependent on the error 


signal and II is dependent on the current through LI. The first three elements 
model the term while the fourth element accounts for the unity term. The 
rest of the model poses no new problems. 





Fig. 7 .6 Controller Circuit Model 


To "be specific then we are given the transfer function of the plant 


where 


0(b) = 




-sT 


{\s + l)(T 2 s + 1 )(T 3 s + 1) 


a 1.0 

T x a 13.1 min 
Tg = 11.1 min 
Tg = 0.5 min 
T a 9*5 min 


(7.18) 


To obtain the unit step response we must determine the controller gain K c and 
the reset time T^. From [LU 1] we have the experimental results obtained with 
an analog computer. This is shown in Fig. 7.7 where the key response charac- 
teristics are labeled and the controller parameters used by Lupfer and Oglesby 
are indicated. 

Following- Eiseriberg [EX 2] we can restate the design problem as: 

Choose values of controller gain constant K e and reset 
time T^ to yield an output time domain response with 


P?1 





Fig. 7.7 Reactor process; response to unit 

step input, Kq - 1.0 and T r - 20.8 min. 

the characteristics of Fig. 7 . 7 . Specifically, it is desired 
that (l) the peak overshoot = 27 $, ( 2 ) the time to reach a 
peak Tp = ^5 min (where the output response begins after the 
initial T = 9.5 min delay), and (3) the time for the response 
to reach 2 $ of the steady- state value T s = lkO min. 

At this point we do not have to reformulate the system transfer function as in 
[EX 2] but simply obtain the UASAP model shown in Fig. 7.8. 

Finally we obtain the unit step response using NASAP. Compare this out- 
put Fig. 7*9 with that given in Fig. 7 * 7 * Ihis ability to handle transport 
lag considerably broadens the range of aerospace /1 control problems that can be 
assisted by the use of KASAP. To permit detailed comparisons Fig. 7.9a is based 

on K =1 and T = 20.8 min. while Fig. 7.9b is a rerun of this example using the 
c r 

values K =0.9 and T =21.9 min. taken from [El 2]. 
c r 
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Fig. 7*8 NASAP Model for Fig. 7-5 
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Vila . BXAMHiE ' SHOWING NASAP LIMITATION - 

We now present an active filter circuit wmcn illustrates a xiowgrapn 
of great complexity and demonstrates some NASAP limitations. In particular 
the fact that it may not be sufficient that the user be judicious in the 
specific tree he allows the NASAP program to select. In Fig. 7-9 
which is identical to Fig. 2.16 in [HU 1] is shown the 28-element NASAP 
circuit diagram of this Hutton problem. The input voltage Vl is fed through 
a low-pass tee-network to the base of a transistor which is part of a two- 
transistor differential amplifier. The output voltage of this differential 
amplifier is connected directly to the base of a simple common emitter 
transistor stage. The voltage at the collector of this transistor, which is 
also the output voltage of the circuit, is fed back through, an RC twin-tee 
network to the base of the other transistor in the differential amplifier. 

Each of these three transistors is represented in Fig. 7-9 b by the h-parameter 
equivalent circuit with h^" 0 and with a capacitor (C^y C^, C^) connected 
between the base and collector terminals. This capacitance is included to take 
into account the frequency characteristic of the transistor . 

The NASAP input listing used by Hutton is reproduced as Fig. 7*10. 

There aid;-’ 17 - of the resistors are listed first and followed* by the seven 

capacitors, in numerical sequence. With this listing, the NASAP tree selection 

^ \ - 

algorithm selects as branches of the tr»ee..\ the elements in the following order; 
VI, R17, Cl, C2, C3, C4, C6, C7, R4, R5, R8, R15. This particular tree generates 
a flowgraph possessing a total number of loops of all orders of 2,440,105 
a very complicated flowgraph indeed. As noted 1 by Hutton, 10 minutes of ex- 
ecution time were required on the UNIVAC 1108. 

Utilizing one of the options discussed in Chapter II, a tree can be 
selected to yield a flowgraph with considerably fewer loops. Seven of the 
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Fig. 7.9^b Hutton Problem 
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elements (namely, VI, R17, Cl, C3, C6, R8, R15) in the above described tree 
are connected to the node numbered 1 in Fig. 7*11 . However, three resistors, 
R6, R7, and R13 are also connected to node 1. If these three elements can be 
made branches of a tree which also contains the seven elements from the tree 
selected by Hutton, then ten of the twelve necessary tree branches will be 
connected to a common node. Such a tree with a definite star-like structure 
should'yield a flowgraph with fewer loops. We shall now indicate how R6, R7, 
and R13 can be selected as tree branches and which branches of the original 
Hutton tree must be removed to make way for these resistors. 

The resistor R7 forms a closed path with the original tree branches C7 
and R17. Since the voltage across R17 is the specified output variable, R7 
will become a tree branch only if C7 can be removed from the tree. This 
is easily accomplished by adding a dummy voltage source dependent on the 
voltage across R7. Likewise the resistor R13 forms a closed path with the 
original tree branches C3 and C4. Since C3 is connected to node 1, we wish 
to make R13 a tree branch in place of C4. A dummy voltage source dependent 
upon the voltage across R13 easily accomplishes this by making R13 a type 2 
element instead of a type 4 element (see description of tree selection 
algorithm) . Finally R6 can be included in the list of tree branches by 
making it the first resistor described in the NASAP input listing. With R6 
in the tree, either R4 or R5 must be removed from the tree. The choice is 
easily made by noting that the co-tree element R3 will form a closed path 
with R4 and C2 when R4 is a tree branch. On the other hand, with R5 
as a tree branch, R3 will form a closed path with R6, R7, R5 and C2. Hence 

R4 should, remain as a tree branch. This is achieved by having R4 precede R5 

r* * r 

in the NASAP input list. 

The revised NASAP input listing that yields a tree with 10 elements 
connected to the same node is shown in Fig. 7*12. The controlled sources V2 
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•Fig. 7.11 Alternative NASAP Model for Hutton 
(compare Fig. 7.9) 
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and. V3 are the dummy voltage sources .necessary to include R7 and R13 
in the tree selected by NASAP. With these 2 elements included, there are 
30 elements in the equivalent circuit (NOTE: 30 elements are the maximum 

number of elements that can be used on the RCA Spectra 70 and IBM 360 machines 
because of their 32 bit computer words). Note in Fig. 7*l£ that the flowgraph 
generated by this new tree contains 415, 515 loops of all orders. Thus a 
saving of over 2 million loops has been achieved by careful selection of the 
tree. The subsequent execution time was 25 minutes on the Spectra 70/46 
(equivalent to approximately 3 minutes on the IBM 360/75). Note also the 
difference between the transfer function found by use of the original Hutton 
tree (Fig. 7VI0) and that found by use of the tree described here (Fig. 7*12), 
The extra 2 million loops results in considerable error in the coefficients 
of the transfer function, see the discussion in [SE 1]. Note in Fig. 7«12' 
that only the zeros of the transfer function are given. Due to excessive 
floating point overflow, the MULLER (root-finding) subroutine was unable to 
determine the poles of the transfer function. 

Through the cooperation of Prof . Alan B. Macnee the Hutton problem was 
run using : CIRAN (a -program based' on. sta.te variables) on the University of 
Michigan IBM 360/67 computer. The impedances were scaled by 10 , ^ and 
frequency 10 The circuit was found to exhibit a pair of dominant conjugate 
complex poles near 2.4 KHz. For a + 10% change in the value of R13- the Q 
of this pole ranged from 7.5 to 20. The fact that these three analyses took 
only 12.13 seconds of CPU time indicates a severe limitation of NASAP for this 
class of problems. 
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In Fig. 7*13 the compensation is seen to consist of two identical RC lag net- 
work. The transfer function of the RC lag network is 


S 



or 


G 

S + G 


The parameter G, the RC lag network pole, is chosen different from all 
plant poles, usually farther inside the left half plane than any pole of the 
desired system. 



Fig. 7.13 Compensated Third Order Control System 


Fig. is equivalent to Fig. 7*13 with respect to the transfer function 



Fig. 7.14 Equivalent Feedback Control System Representation for Fig. 7.13 





The closed-loop transfer function for y/x shorn in Fig . 7.14 is 


a(S+G) 2 A 3 

y f (S £ +SGX2*-)% 2 (l4¥^'f)' (S 3 +AgS ^A^S+Aq ) 

x = ".aA (fSii+G(2f+e;s+G^(d+e+f; ; (7.21) 

1 v + s i .7 

(S 2 +SG(2+1j )+G 2 (1+B>c ) ) (S 3 +A 2 S 2 +A 1 S+A 0 ) " 

Equation (7 .20) for the desired system, when multiplied top and bottom by 

* - '2 

(S+G) becomes 


y B 3 (S+G)2 

X S 5 +(B 2 +2G)S 4 +(B 1 +2GB 2 +G 2 )S 3 +(B' 0 +2GB 1 +G 2 B 2 )S 2 +(g 2 B 1 +2GB 0 )S+G 2 Bq 

(7.22) 

This equation can be compared :ter^is:e_ with (7*21^ after ; ijLhas been .xa* ?r i-tte ; . 
to clear fractions. Thus we - ; find that aA^ = and (Successively 

(B 2 - A 2 } - 

b = — where G is positive. 


- V + (6 - Ag) (Bg - Ag) 
~ ^ ~ 


with G positive . 


( B 0 - A 0 ) + (B x - A x ) (2G - Av;) + (Bg - A g )( (G-A g f - A^ 
:B 3 



(V A l)(-3G 2 -A 1 +2q^9‘ (Bg-Ag) (-2G(G-A 2 ) 2 + AjAg-A^ 
— 


and finally 


(B 1 -A 1 ) (G 3 -G 2 A g +GA 1 -A 0 ) + (B g -A g ) (G^+A^-GAQ+G^C^Agj^GA^Ag) 


2,„ . x2 


d = 


„ 2 _ 

G G„ 
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VHD EXAMHiE IMPLYING LUEMBERCER OBSERVER 


For the last problem we consider an application of modern control theory 
based on state variables to show the versatility of MSAP. In particular we 
use the Luehberger observer method to implement compensation for a control 
system situation in which all the states of the system are not measurable. 

For present purposes it is sufficient to note that the Luehberger observer 
[LU 1] is a device which constructs an estimate of the system state vector based 
upon the available system inputs and outputs. Then, based upon the reconstructed 
state vector, simple matrix algebra manipulations can be utilized to obtain es- 
timates of the missing states or a combination of the missing states. 

Luehberger has shown that for an n-th order system with m measureable states, 
the order of the required observed need only contain n-m poles-. 7 Furthermore 
these pole locations are arbitrary as long as they are different from the eigen- 
values of the system matrix. 

She Luehberger observer accomplishes the desired result by adding dynamics 
in the feedback path of the control system. 

This theory suggests a unique form for general compensation of third order 
systems wherein the designer can place the closed-loop poles at any desired lo- 
cation. Shis development is adapted from Newman [HE 2]. 

The open loop descriptions of the control system shorn in Fig. 7.13 is 


y 0 

u 


S 3 + A 2 S 2 + AjS + Aq 


( 7 . 19 ) 


For convenience the desired closed-loop system transfer function is expressed as 


y 

x 


S 3 + B 2 S 2 + B-jS + B Q 


( 7 . 20 ) 
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These six equations define unique values for the symbols a, b, c, f, e, and 


d respectively. 

It is worth noting [ME 2] that this method of general system compensation 
may be extended to higher order systems . An N-th order system will require a 
string of N-l RC lag networks. Each lag network is driven by the signal e and 
the signal y, both passing through gain blocks as in Fig. 7.13. The z signal 
is formed by adding the output of the string and a signal equal to fy, as in 
Fig. 7. 13 * The a block is used in cascade with the plant as in Fig. 7.13 . 

It is found easier to design compensators for plant transfer functions 
with numerator polynomials by reducing the coefficients in the equations to 
numbers, instead of trying to derive the general relationship. 

For an illustrative example, we consider the open loop and desired closed 
loop transfer functions for the control system in Fig. 7.13 . 

Open Loop System 

A3 ^ 

^ ' S 3 +A 2 S 2 +A 1 S+A 0 S(S+1)(S+10) (7.24) 


Closed Loop System (Desired) 


z = ?3 j / JL 

x S 3 +B 2 S 2 +B 1 S+B 0 ' s+5 


_20r_ 

S+20\ 


(7.25) 


Following the procedure outlined above, we 
c, d, e, f, as listed. 


determine the constants a, b. 


a = 200 

b = 2.24 


c = 1.7 

d = 4.256 

(7.26) 

e = -14.018 

f = 10.762 
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The parameter G, was chosen as twenty (G = 20) which is different from all 
plant poles. 

Having thus specified the Lueriberger observer we need the corresponding elec- 
tric circuit model for the compensated third order control system. This model 
is shown in Pig. 7» 15 • A HASAP run was made for the step response of this model. 
We note that the transfer function shows six critical frequencies near s = -20; 
two zeros and four poles* The step response shows a slight steady state error. 

Execution time on the RCA Spectra 70/46 was 42.86 seconds. A comparison 
run was made of this third-order compensator system using CSMP (Continuous Sys- 
tem Modeling Program) on the IBM 360/75. The CSMP step response for this 
control system checked very closely except that the execution time was 32 seconds 
(22 seconds CPU). This represents a significant cost advantage for HASAP since 
the 360/75 is faster by a factor of approximately 8 over the Spectra 70/46. The 
NASAP printout is shown in Fig. 7.16. 




Fig. 7.15 NASAP Model for Third Order Control. System in Fig. 7.13 
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APPENDIX A 


STEP AM) RAMP RESPONSE CAPABILITY FOR NASAP 


The algorithm used by MSAP to determine the impulse response of a 
circuit can be easily extended to find the step and- ramp response. Listed 
-“bellow- are "modifications that are necessary to accomplish this. Note 
'"that if, for example, both the impulse and step responses are required for 
a particular circuit, the problem must be executed twice . 

To the COMMON/FREQ/ cards of subroutines MSINP and BODE must be added: 

,KRESP 

The following additions must be included in subroutine MSINP. 

After card 1010, three cards must be added: 


Go to 3 

824 WRITE (6, 117) 1011 \ 

117 FORMAT (35H INCLUDE TYPE OF TIME RESPONSE CARD) ±01'2 

1013 

After card ll80, nine cards must be added; 

READ (5,101) INLST ll8l 

WRITE (6,107) INLST 1182 

KRESP = 1000 , . 1183 

IF (XCHAR (INLST (l ) ) . EQ. 49 . AND TCHAR(lNLST(2) )EQ.54,f KRESP' = 1'. 1184 

IF (ICHAR (INLST (l) ) . EQ. 55- AND j i(tHAR(INLST( 2) ) .EQ.56)KKESP/=s -l~-~ 1185 

IF (iCHAR (INLST (1)). EQ. 62) KRESP = 2 1186 

IF (ICHAR (INLST(I)). EQ. 59 AND ICHAR (INLST (2) ) .EQ.4l KRESP = 3 1187 

IF (KRESP. EQ.1000) GO TO 824 II89 

WRITE (6,ll6) ' ' 1189A 


The following changes must be made in INV: 

Change card 14930 

From DIMENSION F (100), T(lOO), COE (2) 

To DIMENSION F (100), T(lOO), C0E(2), TYIE (2,3) 

Change card 14990 

From WRITE (6,50) . 

^Tp " WRITE (6,50) (TYPE (j, KRESP), 2) 

Change card 15 00 " 

From 50 FORMAT (26HI IMPULSE RESPONSE FUNCTION //7H F(t)=) 

To 50 FORMAT (lEL, 2A4, 18 h RESPONSE FUNCTION //7H F(t)=:) 

Change card 15280 

From WRITE (6,40) TRNS, (T(K), F(k), K=l, 51) 

To WRITE (6,40) (TYPE(J,KRESP), J=l,2), TRNS, (T(K), F(K), K=l, 51) 
Change card 15290 

From 4Q F0RMAT(17H1 IMPULSE RESPONSE // 5H TIME, l8x, 10Al/(E12,4,E23 .8) ) 
To 40 F0RI4AT(lHl , 2A4, 9H RESPONSE // 5H TIME, l8x, 10A1/-GE12.4,E23 .8)') 
Change ccar.dri5300V^ 

From %fcITE' ’(€>41) ' 
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To WRITE (6,41) (TYPE(J, KRESP) , J=l,2) 

Change card 15310 

From 41 FORMAT (iHl, 45X, 17HIMPULSE RESPONSE, /) 

To 41 FORMAT (IHl, 45X; 2A4, <?H RESPONSE, /) 

The following additions must be included in subroutine AINIV. 


After card 14930, two cards must be added: 

DATATYPE /4HIMPU,4HLSE, 4HSTEP, 4H 4H 14931 

C0MM0N/FREQ/THI, FEQ (3), CC,q, KRESP 

After card 15060, twelve cards must be added: 

Go To (77,76,78,76), KRESP 15061 

76 IAD2=IAD2+1 15062 

RootR ( IAD2 , 2 ) = 0 15063 

Rootl (IAD2, 2) =0 15064 

Go To 77 15066 

73 IAD2-IAD2+1 15067 

RootR(lAD2,2) = 0 15068 

Rootl(lAD2,2) = 0.001 15069 

IAD2= IAD2+1 15069 A 

RootR ( IAD2 , 2 ) - 0 15069B 

RootX (IAD2,2) = -0.001 15079C 

77 CONTINUE 15069E 


When the above modifications are incorporated in the NASAP program 
package, the only change from a user's point of view is that he must 
include one card, immediately following the NASAP PROBLEM card, on which is 
punched (beginning in column one) one of the four following comments: 

NO RESPONSE 
IMPULSE RESPONSE 
STEP RESPONSE 
RAMP RESPONSE 

The choice of which comment is to be used depends on what type of time 
response, if any, is desired. ' In actuality, only the first, two letters of 
each of the above comments are really necessary . In cards 1181 through 1189A 
of .subroutine NASINP, the program reads the card after the NASAP PROBLEM card 
and then prints it. If the letters in columns one and two are on R and a A 
respectively, the time response variable KRESP is set equal to 3. If there is 
an N in column one and an 0 in , column two, KRESP thenis set equal to 1. If 
there is an I in column 1 and an M in column two KRESP is set equal to 1; If 
there is an S in column one, KRESP then has a value of 2. If the program 
encounters none of the above letter combinations on the second input card, 
it will print out the error message given on card 1013 and 
about the job. 

The variable KRESP is used in subroutine INV to select the proper elements 
of the two-dimensional matrix TYPE such that the headings of the printout of 
the time response agree with the information on the second input card. Thus., 
if a step response is desired (KRESP—2 ) , the elements in the second column .of 
TYPE will be printed. The variable KRESP is also used in INV to select the pro- 
per poles that must be included with the transfer function poles to yield the 

particular time response. On lines 15061, if KRESPfl( impulse response),- the 

£ %0 / 
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APPENDIX B 


Suggested Revision that combines- Subroutines SENS & SENSS 
of NASA? 69/1 into a single shorter Subroutine SENSS 


In the University of Pennsylvania monthly report for July 1969, 
the possibility was mentioned of eliminating subroutine SENS by using 
the results of the calculations of subroutine SENSS (namely S^.) and the 
four equations. 




■ H-i ' ““S 


H 

= Re s x - 

/Iran , 
VReH 

■*«s: 

s x ' 

■ Re s x + ( 

'Reg . 
v. ImE 

* ImS?^ 
a y 

,H . . 

*x + 3 

H 

ImS x 



1 « 3fl 

= Re H + 3 

Im H. 



where 


This can be accomplished by malting the following additions to 
subroutine SENSS. After card 17000 , add the following 3 cards: 
cards 15370, 15380, and 15390 from SENS. After card 17030, remove 
cards 17040, 17050, and 17060 and add the following 6 cards: 

COMMON/LOOPS /DM( 100 ) , PHIH(250), LOC-F(250), LGBSNS(250), ABSH(250), 

1 ABSENS(250), REH(250), IMHi(250), HHSNS(250), CSENS(250), SENSRE(250), 


card 15440 from SENS 

card 15450 from SENS with DU!£L(288) changed to DUMl(38) 
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card 15^90 from SENS 

card, I 550 O from SENS -with the addition of , LGBSNS . 

After card 17110 , add the following 2 cards: cards 15520 and 
15530 from SENS . 

Ai'ter card 17540, add the following 39 cards: 
card 15540 from SENS 

card 15550 from SENS with 1H0 changed to 1H1. 
card 15560 from SENS 
Do .802 1=1, NITR 

Z = 10.** LOGF(t) 

BE = BEAL (CSENS(I)) 

AIM = AIMAG (CSENS(l) ) 

CALL QZERO /REH, I , SENSBE , I , GLAC-, $ 8 l, $ 82 ) 

81 SENSRE .(I) = HE - (IMH(I)/REH(I) ) * AIM 

82 CALL Q£SR 0 ( IMH , I , SENSIM, I , GLAG , $ 83 , $84) 

83 SENSIM (I) = BE -r (REH(l)/lMH(l)) * AIM 
■84 SEN 8 BS (I) = RE- 

97 CALL QZERO ( PiU-H, I , SENSEI , I , GLAG, $ 98 , $ 99 ,) 

98 SENSFI (I) = (180. /PI) * (AlM/PHIH(l) ) 

99 WRITE (6,180) L0GF( I ) , Z , SENSBE ( I ) , SENSIM( I ) , SENSBS ( I ) , SENSFI ( I ) 
180 FORINT ( 5 X, 6 (E 15 . 7 , 2X) ) 

802 CONTINUE 

cards l6ll0 thru 16320 from SENS (22 cards) 
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The following 10 cards of SENSS must he slightly modified: 
(NOTE - etc. means that the rest of the card remains unchanged). 
Change card 17120 
from NRITE(6,100) 
to WRITS(6,200) 
card 17130 

from 100 F0RMAT(1H1, etc. 
to 200 FORMA T( 1H0, etc. 
card 17390 

from 9 6 WKITE(6.,105) etc. 
to 9 6 WRITE (6,205) etc. 
card 17400 

from 105 FORMAT etc. 
to 205 FORMAT etc. 
card 17450 

from WRITS (6, 102) 
to WRITE (6,202) 
card 17460 

from 102 FORMAT etc. 
to 202 FORMAT etc. 
card 17470 

from WRITE(6,103) etc. 

;o >©112(6,203) etc. , 

{-25k. y ' 



card 17480 

from 103 FORMAT etc. ' 
to 203 FORMAT etc. 
card 17520 

from WRITE(6,102) 
to WRITE (6,202 } 
card 17530 

from WRITE(6,103) etc. 
to WRITE(b,203) etc. 

Finally in BODE remove card 13510. In SPLOT the following 2 minor 
modifications are needed: 

On card 16670 

change ARR(lOOO), DUMPS (17 50) 
to ARE (750) , DUMM2 (2250 ) 

On card 16690 

change DUML(288) 
to DUML( 38) 

These changes essentially delete cards 15570 thru l6l00 from SENS. 

The calculations in this block of cards are replaced by the Do-Loop 
(labelled 802) given above. The remaining cards of SENS are then inserted 
into subroutine SENSS in the appropriate locations. On the RCA Spectra 70, 
the original version of SENS requires 4000 bytes of memory and fche origi- 
nal version of SENSS requires 2540 bytes. However 'Che suggested new 
version of SENSS requires 4392 bytes - thus realizing a saving of 21^8 
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bytes. Enclosed is a listing of the new version of SENSS . Due to the 
use of a BCD coded program on a EBCDIC machine, the following symbols 
are equivalent: 

^ is equivalent to = 

% " " " ( 

< M n " ) 

n »» « » . 

CS - V 

Note that the printed output of the original SENSS now precedes 
; the printed output of the original version of SENS. 
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